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influence?  A  and  thu?  affect?  the  prediction?  of  A  .  U?ing  critical  ?tate 
concept?  a  simplified  procedure  to  predict  the  <-Ver ocnsciidatio  ratio  of  clay? 
i?  presented.  Since  the  OCR  estimate  i?  not  ?en?itivc  to  the  C  /C  ratir,  a 
good  prediction  i?  po??ible  from  the  knowledge  of  e/e  and  oh  c 

The  limitation?  of  existing  theories  to  model  the  ani?ctropy  of  clay?  were  .iden¬ 
tified.  A  novel  approach  I?  presented  that  incorporate?  the  essential  feature? 
of  anisotropic  behavior  of  clays,  and  the  theory  is  validated  using  experiments* 
data.  The  model  capture?  the  essential  features  of  initial  and  induced  aniso¬ 
tropy,  yet  yielding  close  form  solution?  for  all  the  failure  parameter?  of  a  ay. 
Tr.e  propoped  theory  i?  applicable  to  a  number  of  geotechnical  problem?;  ir.  th*s 
report,  it  i?  used  to  investigate  several  important  issues  related  to  the 
determination  of  in  situ  properties  using  the  seif-boring  pressuremeter  ( S5PK : : 

(i)  The  stress  state  and  mode  of  deformation  ir.  pressuremeter  tests  are- 
identified.  The  generality  of  plane  strain  condition  during  expan¬ 
sion  of  an  SBFM  cavity  is  shown  ar.d  condition?  where  it  is  far  from 
plane  strain  are  indicated.  Experimental  data  are  used  to  substan¬ 
tiate  the  tnecre tica*  findings.  Ail  indication?  are  that  these 
finding?  car.  settle  the  questions  regarding  the  mode  of  deformation 
and  the  role  of  the  vertical  stress  in  expansion  of  cylindrical 
cavity. 

(ii)  Tnc  possibility  of  radial  cracking  (negative  circumferential  stress) 
war  analyzed  ir  detail  and  shown  to  be  unlikely. 

(1*1  )  U?i.ng  the  theory  developed  ir:  this  work,  the  SB  PM  failure  parameters 
car  be  used  to  determine  the  same  parameter?  for  many  different  mode? 
of  failure  other  than  that  of  the  SBPK. 

(ivi  A  procedure  ha?  beer:  proposed  to  evaluate  the  error  in  strain  induced 
by  initial  movement  of  the  pressuremeter  cavity.  After  evaluation, 
of  the  error  ir  strain,  actual  strength  parameters  car  be  calculated. 
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CHAPTER  1 
INTRODUCTION 

The  research  study  undertaken  at  Purdue  focuses  towards 
increasing  our  understanding  of  fundamental  issues  related  to  the 
behavior  of  clays,  especially  with  regard  to  modelling  their  aniso¬ 
tropy.  The  information  being  developed  is  applicable  to  most 
geotechnical  problems,  however,  in  illustrating  their  usefulness, 
emphasis  is  placed  on  interpretation  of  in  situ  tests,  and  in  par¬ 
ticular  pressuremeter  and  self-boring  pressuremeter  (SBPM)  tests. 
There  is  ample  evidence  that  in  situ  soils  are  anisotropic, 
e 1  a s t o-pl a s t i c  ,  stress  path  and  rate  dependent,  unlike  what  is  gen¬ 
erally  assumed  in  te6t  interpretations.  Therefore,  the  study  of 
mechanisms  of  deformation  of  in  situ  tests  and  their  interpretation 
must  take  into  account  these  features,  and  especially  the  anisotro¬ 
pic  nature  of  the  soil  deposit.  A  simple  and  reliable  anisotropic 
theory  will  be  most  useful  in  the  study  of  several  important  fac¬ 
tors  related  to  pressuremeter  testing  such  as  possibility  of  radial 
cracking,  role  of  vertical  stress,  stress  conditions  at  failure, 
and  the  effects  of  initial  disturbance  which  cannot  be  measured. 

It  will  also  be  possible  to  estimate  undrained  strength  of  the  clay 
in  other  modes  of  failure  using  SBPM  data  only. 

In  Chapter  2  of  this  report,  the  viability  of  modelling  the 
anisotropic  behavior  of  clays  within  the  context  of  the  modified 
cam  clay  (MCC)  model  (Roscoe  and  Burland  1968,  Wroth  1984)  is  stu- 


died.  Some  useful  relationships  to  predict  the  behavior  of  clays 
in  triaxial  space  are  presented,  and  the  limitations  of  the  MCC 
model  are  indicated.  A  cuboidal  shear  device  is  used  to  illustrate 
the  influence  of  stress  path  and  fabric  anisotropy  on  the  behavior 
of  clays.  In  addition,  using  critical  state  concepts  (Schofield 
and  Wroth  1968,  Wroth  1979)  a  simplified  procedure  to  predict  the 
o ve r c on s ol i da t i on  ratio  of  clays  is  presented. 

In  Chapter  3,  the  limitations  of  existing  theories  to  model 
the  anisotropy  of  clays  are  identified.  A  novel  approach  is 
presented  that  incorporates  the  essential  features  of  anisotropic 
behavior  of  clays,  and  the  theory  is  validated  using  experimental 
data.  The  model  captures  the  essential  features  of  initial  and 
induced  anisotropy,  yet  yielding  close  form  solutions  for  all  the 
failure  parameters  of  clay.  All  the  facets  of  anisotropy  known  to 
date  are  incorporated.  It  can  be  applied  to  a  number  of  geotechni¬ 
cal  problems,  such  as  stability  of  embankments,  retaining  systems, 
etc.  Its  potential  is  demonstrated  by  studying  the  following 
issues  in  p r e s s u r e me t e r  testing  of  soils.  (i)  The  stress  state  and 
mode  of  deformation  in  p r e s s u r e me t e r  tests  are  identified.  The 
generality  of  plane  strain  condition  during  expansion  of  an  SBPM 
cavity  i6  shown  and  conditions  where  it  is  far  from  plane  strain 
are  indicated.  Experimental  data  are  used  to  substantiate  the 
theoretical  findings.  All  indications  are  that  these  findings  can 
settle  the  questions  regarding  the  mode  of  deformation  and  the  role 
of  the  vertical  stress  in  expansion  of  cylindrical  cavity  (SBPM). 


(ii)  Possibility  of  radial  cracking  is  analyzed  in  detail,  showing 
little  chance  for  its  occurrence.  (iii)  A  method  to  obtain  aniso¬ 
tropic  failure  parameters  using  SBPM  test  data  is  outlined.  (iv) 
An  analytical  technique  to  make  allowance  for  the  disturbance  in 
te6t  data  is  presented. 


CHAPTER  2 


STRESS  PATH  AND  STRENGTH  ANISOTROPY 
2.1  Experimental  Illustration  of  Anisotropy  Effects 

Strength  anisotropy  consists  of  two  major  components.  One  is 
inherent  anisotropy,  which  occurs  due  to  preferred  particle 
arrangement  during  sedimentation,  resulting  in  a  fabric  that  is  not 
identical  in  all  directions.  The  other  component  is  stress  induced 
anisotropy  which  is  caused  by  an  anisotropic  state  of  stress  at  the 
end  of  consolidation.  Inherent  anisotropy  implies  that  the  soil 
behaves  anisotropically  even  if  the  initial  stress  state  is  isotro¬ 
pic.  Stress  induced  anisotropy  means  that  the  6oil  behaves  aniso¬ 
tropically  depending  on  the  direction  of  loading  due  to  its  initial 
anisotropic  stress  state  even  if  the  soil  properties  such  as  c'  and 
<t>'  are  isotropic.  The  ratio  of  triaxial  compression  strength  to 
that  of  triaxial  extension  is  often  referred  to  as  a  measure  of 
anisotropy.  Using  the  cuboidal  shear  device,  anisotropy  was  illus¬ 
trated  by  comparing  the  strength  derived  from  strain  controlled 
tests  on  sedimented  samples  loaded  vertically  and  horizontally 
under  plane  strain  conditions. 

Three  tests  were  carried  out  on  WP  Georgia  silty  clay  under 
essentially  plane  strain  conditions:  the  first  two  tests  on  Kq 

consolidated  samples  and  the  following  one  on  an  isotropically  con¬ 
solidated  sample.  In  the  first  test  (CKoUC-PSV),  major  and  minor 
principal  stress  increments  were  applied  horizontally,  while  the 
Intermediate  principal  stress  increment  was  applied  vertically. 


ensuring  plane  strain  conditions,  with  the  plane  of  deformation 

being  horizontal  (this  situation  is  typical  of  a  horizonta  Iiy 

loaded  pile  foundation).  In  the  second  test  (CK  UC-PSH)  intermedi- 

o 

ate  and  minor  principal  stress  increments  were  applied  horizon¬ 
tally,  while  the  major  principal  stress  increment  was  applied  vert¬ 
ically,  ensuring  plane  strain  conditions,  with  the  plane  of  defor¬ 
mation  being  vertical  (this  is  a  typical  field  loading  condition, 
e.g.,  long  embankment  strip  footing,  etc.).  The  same  test  was 
repeated  on  an  isotropically  consolidated  sample  (CIUC-PSH).  The 
results  of  these  plane  strain  tests  are  given  in  Table  2.1. 

To  quantify  strength  anisotropy,  several  techniques  have  been 
proposed  by  previous  researchers  (e.g.  Aas,  1965,  Duncan  and  Seed, 

1  966  ,  Lo  and  Morin,  1  9  7  2  ,  Berre  and  Bjerrum,  1  9  7  3  ,  K r i s h na mu r t h y  , 
et  al.,  1980,  Nakase  and  Kamei,  1983).  In  the  experimental  program 
described  herein,  the  specimen  is  never  rotated.  The  axis  remains 
vertical  during  consolidation  and  shear.  Anisotropy  is  measured  by 
the  ratio  of  the  shear  strength  of  a  horizontally  loaded  specimen 
to  the  strength  of  a  vertically  loaded  replicate  specimen. 

From  Table  2.1,  the  three  different  tests  that  followed  dif¬ 
ferent  stress  paths  and  had  undergone  different  consolidation  (K^ 
and  i6ot  ropic  consolidation)  yielded  different  strength  values. 

The  first  sample  (CK  UC-PSV),  K  consolidated  and  loaded  horizon- 

o  o 

tally,  gave  the  lowest  strength.  The  reason  is  that  the  extension 
test  Btress  path  meets  the  failure  envelope  before  the  compression 
test  6tress  path  does.  In  the  other  two  tests  where  the  plane 


Table  2.1 

Plane  Strain 

Test  Result 

(WP  Georgia 

Silty  Clay) 

Sample 

Test 

(Tf/0cv)max 

1  V°cv>H 

* 

No. 

1 

CK  UC-PSV 

0.47 

0 

0.78 

2 

CK  UC-PSH 

0 

0.  60 

3 


CIUC-PSH 


0.  70 


strain  condition  was  maintained  in  horizontal  direction  and  the 

loading  was  applied  vertically  (conventional  plane  strain  tests), 

the  samples  showed  higher  strengths  th3n  the  first  one,  with  the 

isotropically  consolidated  sample  being  stronger  than  the  K  conso- 

o 

lidated  one. 

The  three  different  stress  path  tests  yielded  different 

strengths  which  clearly  illustrate  the  6tress  dependent  anisotropic 

nature  of  the  soil.  The  anisotropic  ratio  [t,/o  ]„/[t,/o  ] „  of 

fvcH  fvcV 

the  Kq  consolidated  sample  is  0.78,  which  shows  that  the  stress 
path  has  a  significant  influence  on  the  anisotropic  strength.  The 
strength  ratio  between  the  Kq  and  isotropically  consolidated  sam¬ 
ples  is  not  unity  ( ■  0.86)  which  indicates  that  induced  anisotropy 
effects  exist. 

2.2  Theoretical  Studies  of  Stress  Path  and  Anisotropy 

2.2.1  Predicting  Undrained  Shear  Strengths  of  CK  UC  from 
CIUC  Tests  ° 

Isotropically  consolidated  undrained  compression  (CIUC)  triax- 
ial  tests  are  the  common  laboratory  tests  to  determine  the 
undrained  shear  strength  of  clays  primarily  because  of  the  conveni¬ 
ence  and  simplicity  of  the  experimental  procedures.  Anisotropi- 
cally  consolidated  (CAUC)  tests  are  more  complicated  and  take  much 
longer  to  run.  Natural  soil  deposits,  however,  are  rarely  isotro¬ 
pic  nor  do  they  often  have  an  in  6itu  Kq  •  1.  Therefore,  develop¬ 
ment  of  a  theoretical  expression  to  predict  the  actual 
behavior/strength  (i.e.  CKoUC  strength)  from  CIUC  test  results 


would  be  very  beneficial 


Although  Ladd  (1965)  concluded  that  it  is  not  possible  to 
predict  CK^UC  strength  from  C1UC  tests  because  the  test6  follow  two 
different  stress  paths,  Sivakugan ,  et  al.  (  1  988b)  found  that  a  rea¬ 
sonably  good  estimate  of  the  CK^UC  strengths  may  be  obtained  from  a 

single  CIUC  test  on  the  same  soil.  The  ratio  of  (t^/g^  ) CK  uc  t0 

o 

(t,/o  )„T,.„  was  shown  to  be  a  function  of  the  coefficient  of  earth 

f  vo  CIUC 

pressure  at  rest  and  Skempton's  pore  pressure  parameters  at 

failure  Ar  for  both  tests.  Thus,  T_/o  from  a  CIUC  test  may  be 
f  f  v  o 

used  with  estimates  of  K  ^  and  the  A  parameters  to  estimate  the  nor¬ 
malized  shear  strength  of  K^-consolidated  specimens.  Although  only 
limited  data  are  available,  the  procedure  appears  to  be  suffi¬ 
ciently  accurate  for  practical  purposes.  It  is  briefly  summarized 
here,  as  details  are  given  in  Sivakugan  (1987)  and  Sivakugan,  et 
al.  (1988b). 


It  was  shown  that  for  normally  consolidated  clays,  the 
strength  ratio  can  be  given  by: 


whe  re 


v  'CK  UC 
o  o 

vo _ 

<4i> 


CIUC 


vo 

T 


( - “) 

v  '  'CK  UC 
o  o 

vo 


(-^) 


CIUC 


vo 


K  +2(1  -  K  )  A ,  . 

o _ ' _ Q  f  ,  1 

K  +2(1  -K)Af 

O  O  I  ,  K 

9  A 


Uf  v  (l 
f  *0 


K  )  +  K  ]  (2.1) 

o  o 


normalized  undrained  shear  strength  for  CK^UC  test 
normalized  undrained  shear  strength  for  CIUC  test 
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A,  ■  Skempton's  A  parameter 

t  i  * 

A  -  Skempton's  A  parameter 

1  >  K 


K  -  coefficient  of  at  rest 
o 


at  failure  for  the  CIUC  test; 

at  failure  for  the  CK  UC  test 

o 

earth  pressure. 


Typical  representations  of  this  strength  ratio  equation  are 

given  in  Fig.  2.1  for  Kq  “  0.4.  From  this  figure  and  others  at 

different  values,  it  can  be  seen  that  A^  ^  has  a  much  greater 

influence  on  the  strength  ratio  than  A,  ,  and  even  a  crude  esti- 

f  ’o 

mate  of  A,  is  usually  sufficient  for  a  reasonably  good  estimate 

£*  o 

of  the  strength  ratio. 


The  parameters  required  for  the  computation  are  K  ,  A  and 

o  f  ,  i 

A  .  It  is  admittedly  not  easy  to  obtain  accurate  estimates  of 
*  o 

Kq  in  practice,  e.g.  see  Massarsch  et  al.  (1975)  and  Tava^as  et  al. 
(1975).  In  the  interim,  Kq  may  be  obtained  from  published  correla¬ 
tions  (Holtz  and  Kovacs,  1981)  or  estimated  from  Jaky'6  (1948) 
relationship.  The  isotropic  A  parameter  A  can  be  obtained  from 

*  9  * 

the  CIUC  test  results,  while  Ar  „  requires  K  -triaxial  tests  (lust 

f  ,  K  o  J 

o 

what  we  want  to  avoid).  Fortunately,  as  indicated  previously,  the 

strength  ratio  is  not  very  sensitive  to  A,  ,  and  only  crude  esti - 

*'0 

mates  are  required  for  good  estimates  of  the  strength  ratio.  In 

the  absence  of  any  other  information,  the  mere  knowledge  of  whether 

A,  .  is  greater  or  less  than  A,  v  is  sufficient  for  a  reasonably 
r,i  l,Ko 

good  estimate  of  this  parameter.  Numerical  examples  to  predict  the 
strength  ratio  are  given  in  Sivakugan  et  al.  (1988b). 


To  check  the  applicability  of  the  proposed  technique,  calcu¬ 
lated  values  of  the  strength  ratio  were  compared  to  experimental 
data  (Table  2.2)  for  24  different  soils,  and  the  results  are  plot¬ 
ted  in  Fig.  2.2.  This  figure  shows  that  for  carefully  performed 
tests  on  replicate  specimens,  the  assumptions  made  and  equations 
used  to  derive  Eq .  2.1  are  reasonably  valid.  From  Table  2.2,  it  is 

seen  that  A  can  be  less  or  greater  than  Af  by  as  much  as 

*»  o  f>i 

100%.  The  slight  differences  in  tp'  between  isotropic  and  K  tests 
can  be  neglected  in  the  procedure.  This  is  justified  by  a  statist¬ 
ical  analysis  of  experimental  data  from  C1UC  and  CAUC  tests  by 
Mayne  (  1  985  ),  which  showed  that  4>CAUC  =  0.97  ‘f'  c  j  jj  C  * 

2.2.2  Spacing  Ratio 

The  Cam  clay  (Roscoe  et  al.,  1963)  and  the  modified  Cam  clay 
models  (Roscoe  and  Burland,  1968)  are  the  most  widely  accepted 
critical  state  6oil  models.  They  have  undergone  several  develop¬ 
ments  over  the  past  two  decades  (e.g.,  Egan,  1977,  Pender  1978,  and 
van  Echelen  and  Potts,  1978).  Sivakugan  (1987)  and  Sivakugan  et 
al.  (1988a)  give  a  review  of  the  main  characteristics  of  these 
models.  The  Cambridge  soil  models  were  developed  essentially  for 
isotropically  consolidated  clays,  although  most  soil  deposits 
encountered  in  nature  are  one  dimensionally  consolidated  with  no 
lateral  deformation. 

In  the  Cambridge  soil  models,  it  was  hypothesized  that  the  K 

o 

consolidation  line  (Kq  CL)  in  the  p'-q-v  space  lies  on  the  state 
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Computed  from  Eq .  I  of  Bjerrum  and  Simmons  (I960) 


boundary  surface  between  the  critical  state  line  (CSL)  and  the  iso¬ 
tropic  consolidation  line  (ICL)  (Atkinson  and  Bransby,  1978,  and 
Roscoe  and  Poorooshasb,  1963).  Furthermore,  ICL,  KqCL,  and  CSL  are 
assumed  to  be  three  parallel  lines  with  a  negative  slope  of  A  in 
the  v-ln  p'  plane,  as  shown  in  Fig.  2.3  (Atkinson  and  Bransby, 
1978),  resulting  in  the  following  equations,  respectively: 

v  =  N  -  A  li,  p'  (2.2a) 

v  =  Nq  -  A  In  p'  (2.2b) 

v=r-Alnp'  (2.2c) 

where  N  =  the  specific  volume  (v  =  1  +  e)  of  a  specimen  isotropi- 

2 

cally  (hydrostatically)  consolidated  under  p'  *=  1  kN/m  .  Nq  is  the 

specific  volume  of  a  specimen  one  dimensionally  consolidated  with 

2 

p'  “  1  kN/m  ,  and  T  is  t!  e  specific  volume  at  the  critical  state 

2 

when  p '  =  1  kN/m  . 

A  new  state  parameter,  called  spacing  ratio  and  denoted  by  r, 
is  introduced  herein  (Sivakugan  et  al.  ,  1988a)  to  represent  the 

relative  position  of  KqCL  with  respect  to  ICL  and  CSL.  The  stress 
path  of  a  CKoUC  test  is  AF  (Fig.  2.3),  and  the  spacing  ratio  can  be 
defined  as: 

n  -  r 

r  -  ~---r  ■  (0  <  r  <  1)  (2.3) 

When  this  ratio  takes  its  extreme  values  of  0  or  1 ,  K  CL  coincides 

o 

with  the  CSL  or  ICL,  respectively. 

Expressions  were  derived  for  spacing  ratios  for  extended  Cam 
clay  model  and  for  extended  modified  Cam  clay  model  in  order  to 
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consider  Kq  consolidated  soils  instead  of  isotropically  consoli¬ 
dated  soils.  The  spacing  ratio  for  both  extended  models  are  given 
in  Table  2.3.  Sivakugan  et  al.  (1988a)  have  derived  expressions 
for  the  normalized  6hear  strength  T  /a  using  the  Cam  clay  model, 
modified  Cam  clay  model  and  the  extended  models  with  spacing  ratio. 
The  derived  expressions  are  summarized  in  Table  2.3.  For  detailed 
derivations  the  reader  is  referred  to  Sivakugan  (1987).  The  nor¬ 
malized  shear  strength  is  a  function  of  tne  friction  angle  and  the 
consolidation  characteristics  of  the  soil.  The  expressions  are 
derived  assuming  that  the  normally  consolidated  clay  is  idealized 
as  an  elasto-plastic  material  exhibiting  isotropic  strain  harden¬ 
ing.  Based  on  constitutive  equations  proposed  by  Sekigushi  and 
Ohta  (1977),  Ohta  and  Nishihara  (1985)  developed  similar  equations 
for  the  normalized  shear  strength  but  using  rheological  and  dila- 
tancy  characteristics  of  soils. 

A  typical  variation  of  normalized  undrained  shear  strength 
with  4>'  predicted  by  the  extended  models  is  shown  in  Fig.  2.4.  The 
extended  modified  Cam  clay  model  gives  slightly  higher  values  than 
the  extended  Cam  clay  model  for  all  values  of  Kq  and  A,  however, 
the  differences  are  small.  The  normalized  undrained  shear  strength 
varies  approximately  linearly  with  $  .  It  varies  between  0.25 

L  l  u  w 

and  0.45  for  typical  values  of  Kq  ,  A  and  (Sivakugan  et  al., 
1988a)  . 

For  the  limited  data  available  in  the  literature,  predictions 

of  t , / o  for  K  consolidated  clays  were  made  from  the  extension  of 
f  vo  o 


Table  2.3  Expressions  for  Spacing  Ratio,  Normalized  Shear 
Strength  and  Skempton's  A-Parameter  at  Failure 


Extended  Cam  Clay 
Extandad  Modified  Cam  Clay 


the  Cam  clay  and  modified  Cam  clay  models 


The  data  and  the 


predicted  values  are  given  in  Table  2. A.  Both  extended  models  show 
excellent  agreement  with  the  experimental  results  as  illustrated  in 
Fig.  2.5. 


The  Cam  clay  models  as  well  as  the  extended  models  can  also  be 
used  to  predict  the  pore  pressures  at  failure  by  evaluating 
Skempton's  A  parameter  at  failure,  A  ^  .  This  parameter  is  given  by: 


At  failure  q ^  =  Mp  ^ .  Therefore, 


a  i  i  .  I 

Af  3  MM' 

Pf 


(  2.  A  ) 


(2.5) 


The  resulting  relationships  between  A  _  and  M,  r  r,  and  A, 

to 


are  given  in  Table  2.3  (see  Sivakugan,  1987,  for  a  detailed  deriva¬ 
tion).  A  typical  variation  of  A,  ,  the  A  parameter  at  failure 

t  ,  K. 


for  K  consolidated  clays,  with  predicted  from  the  extended 

o 


models  is  shown  in  Fig.  A.  For  low  values  of  A  both  models 

predict  about  the  same  A  .  For  higher  A,  the  extended  modified 

o 

Cam  clay  model  gives  lower  values  than  the  extended  Cam  clay  model 


A,  decreases  with  increasing  ,  decreasing  A,  and  increasing 
1  ’  o 


K  .  For  typical  values  of  A  and  K  ,  Af  varies  between  0.9 


O  A 


and  4.0  (Sivakugan  eC  al.,  1988a). 


Predictions  of  A,  for  K  consolidated  clays  were  made  for  the 

f  o 

eight  soils  given  in  Table  2.4.  The  data  and  results  are  given  in 
Table  2.4  for  both  the  extended  Cam  clay  model  and  extended  modi¬ 
fied  Cam  clay  models.  Unlike  the  excellent  agreement  obtained  for 

the  T  / o  ratio,  the  A,  prediction  from  both  models  are  incon- 
f  vo  f 

sistent,  being  excellent  for  some  tests  and  poor  for  others  (Fig. 

2.7).  It  seems  that  A^  is  very  sensitive  to  A,  which  is  a  function 

of  both  C  and  C  ,  the  compression  index  and  the  swelling  index, 
c  r 

respectively.  Measurement  of  Cc  is  rather  straight  forward;  how¬ 
ever,  C^  varies  with  OCR  and  no  standard  procedure  is  generally 
adopted  in  measuring  C^.  This  severely  influences  A  and  thus  sig¬ 
nificantly  affects  the  predictions  of  A^.  .  To  avoid  the  problems  in 
determining  A  from  Cf ,  Mayne  (1980)  has  proposed  to  redefine  A  in 
terms  of  o ve r c on s o 1 i da t e d  strength  data  as: 


log  1(Tf/,°v.)0C1  ~  108 
log  OCR 


(2.6) 


where  (t r/a  )_„  and  (t,/o  )„_  are  the  normalized  shear  strength 

f  vo  OC  f  vo  NC 

of  overconsolida ted  and  normally  consolidated  clays,  respectively. 
The  normalized  undrained  strength  values  at  two  different  OCR's  are 
required  for  this  determination.  For  the  data  given  in  Table  2.4, 
such  information  was  not  available,  and  thus  it  was  not  possible  to 
compute  A  by  this  method. 

It  has  also  been  suggested  to  determine  A  from  correlations. 


experimental  A 


O  Extended  Cam  Clay 
•  Extended  Modified  Cam  Clay 


Predicted  Aj 

Figure  2.7  Predicted  vs.  Measured  A-parameter  at  Failure 
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Few  empirical  correlations  have  been  cited  in  the  literature  relat¬ 
ing  M  with  A.  Schofield  and  Wroth  (1968)  proposed  that  M/A  -  1.5; 
Karube  (1975)  suggested  that  it  is  1.75.  However,  based  on  experi¬ 
mental  data  from  literature,  Sivakugan  and  Holtz  (1986)  showed  that 
there  appears  to  be  no  correlation  between  M  and  A.  As  shown  in 
Fig.  2.8,  the  rather  large  scatter  makes  a  linear  correlation  ques- 
t ionable . 

2.2.3  Conclusions 

Although  CKqUC  and  CIUC  tests  follow  different  stress  paths, 

it  has  been  shown  that  a  good  estimate  of  the  CKqUC  strength  can  be 

obtained  from  a  single  CIUC  test.  The  normalized  strength  ratio, 

*  + 

(Tr/°  )  rv  (  T  , /o  )  ,  was  found  to  be  a  function  of  earth 

f  vo  CK  UC  f  vo  CIUC 
o 

pressure  at  rest  Kq  and  Skempton's  pore  pressure  parameters  at 

failure  A^  for  both  tests.  Kq  may  be  obtained  from  published 

correlations  or  estimated  from  Jaky's  relationship.  The  isotropic 

A  parameter  A  can  be  obtained  from  the  CIUC  test  results.  Since 

1,1 

the  strength  ratio  is  not  very  sensitive  to  A  ,  an  estimate  of 

f  ’  o 

A  is  sufficient  for  the  calculation. 

f  ’o 

Expressions  were  derived  for  the  normalized  shear  strength 
t f/°vo  using  the  Cam  clay  model,  the  modified  Cam  clay  model,  and 
the  extended  Cam  clay  and  modified  Cam  clay  models  with  spacing 
ratio.  The  normalized  shear  strength  is  a  function  of  the  friction 
angle  and  the  consolidation  characteristics  of  the  soil.  Both 
extended  models  with  spacing  ratio  predict  values  that  compare  well 


2P 

with  experimental  data. 

Expressions  were  also  derived  for  estimating  Skempton's  A 
parameter  at  failure  for  isotropically  consolidated  soils  and  Kq 
consolidated  soils.  Unlike  the  excellent  agreement  with  experimen¬ 
tal  results  obtained  for  the  strength  ratio,  the  prediction  is 

marginal  for  all  models.  A^  is  very  sensitive  to  A,  which  is  a 

function  of  C  and  C  .  This  severely  influences  A  and  thus  affects 
c  r 

the  predictions  of  A  ^ . 

2.3  Normalization  of  Consolidation  Lines 

2.3.1  Introduction 

A  good  estimate  of  compression  index,  Cc  ,  re c omp re s s i on  index, 
Cf,  and  preconsolidation  pressure,  p^,  are  prerequisites  for  set¬ 
tlement  analysis  of  embankments  or  buildings  on  cohesive  soils. 

For  this,  it  is  common  practice  to  perform  consolidation  tests  on 
specimens  prepared  from  undisturbed  samples  taken  at  various  loca¬ 
tions  within  the  clay  stratum.  From  the  consolidation  test 
results,  e-logo^  plots  are  generated,  and  ,  C^.  and  p^  are 
obtained.  These  steps  involve  effort  and  expense,  and,  for  prelim¬ 
inary  designs,  empirical  correlations  are  often  used  to  estimate 
these  parameters  (e.g.  Holtz  and  Kovacs,  1981).  Cc  has  been  corre¬ 
lated  with  natural  water  content,  liquid  limit  and  initial  void 
ratio  (e.g.  Skempton,  1944,  Nishida,  1936,  Cozzolino,  1961,  Azzouz 
et  al.,  1976,  and  Koppula,  1981).  Cr  is  often  assumed  to  be  0.1  to 
0.2  times  C^.  Estimating  p^  from  correlations  remains  a  problem 
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yet. 


The  slope  and  position  of  virgin  consolidation  lines  in  the 
e-logo^  space  depend  primarily  on  liquid  limit  (or  initial  void 
ratio)  and  sensitivity  of  the  soil.  Therefore,  for  clays  of  low  to 
medium  sensitivity,  the  slope  and  the  intercept  on  the  e-axis 
depend  mainly  on  the  liquid  limit.  When  e-logo^  lines  for  several 
normally  consolidated  clays  are  plotted  together,  a  wide  range  of 
void  ratios  may  be  observed  for  a  given  o^,  depending  upon  their 
respective  liquid  limits.  Based  on  the  Gouy-Chapman  diffuse  double 
layer  theory,  Nagaraj  and  S r i n i va s amu r t hy  (1983  and  1986)  showed 
that  when  e  is  normalized  with  respect  to  the  void  ratio  at  the 
liquid  limit,  e^ ,  all  the  lines  tend  to  fall  into  a  narrow  band 
irrespective  of  the  clay  mineral  type  or  pore  fluid  characteris¬ 
tics. 


This  is  confirmed  from  published  consolidation  test  data  for 
16  different  clays  and  experimental  data  from  three  artificially 
sedimented  clays.  The  existence  of  such  a  normalized  relationship 
facilitates  predictions  of  the  o ve r c on s o 1 i d a t i on  ratio  (OCR).  A 
simple  procedure  to  predict  OCR  is  proposed.  A  summary  of  this 
study  and  it6  conclusions  are  given  below. 

2.3.2  Critical  State  Model 

The  virgin  consolidation  line  is  often  assumed  to  be  linear 

(i.e.  constant  C  )  for  settlement  calculations.  An  idealized  vir- 

c 

gin  consolidation  line  in  e-logo  space  is  shown  in  Fig.  2.9.  In 
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reality,  it  may  be  slightly  curved  as  the  plastic  limit  is 
approached.  The  equation  of  the  virgin  consolidation  line  is  given 


e  *  -  C  logo  +  e  , 
c  6  v  1 


(2.7) 


Normalizing  Eq.  2.7  by  e  -  e  : 

L<  r 


e  "Cc  el 

-  C  -  log  0  +  - 

eL  “  eP  eL  '  6P  V  eL  ‘  eP 


(2.8) 


where  ep  is  the  void  ratio  at  plastic  limit.  The  liquidity  index 
LI  is  defined  bv: 


w  -  PL 
LL  -  PL 


From  Eqs.  2.8  and  2.9 


eL  -  CP 


C  ,  e ,  -  e  _ 

c  ,  _  1  P 

-  log  o  +  - 

p  —  p  V  p  —  p 

L  P  L  P 


(2.9) 


(2.10) 


Thus,  if  e  versus  logOv  is  a  straight  line,  so  is  LI  versus  logOv 
and  vice  versa. 


There  is  experimental  evidence  that  when  the  liquidity  index 

is  plotted  against  logarithm  of  undrained  shear  strength,  all  the 

points  fall  within  a  narrow  band  (e.g.  Skempton  and  Northey,  1953, 

Houston  and  Mitchell,  1969,  Mitchell,  1976).  Based  on  the  large 

number  of  experimental  results  given  by  Skempton  and  Northey  (1953) 

and  Youssef  et  al.  (1965),  Schofield  and  Wroth  (1968)  adopted  a 

2 

value  of  1.7  kN/m  as  the  best  estimate  of  the  undrained  shear 


strength  at  the  liquid  limit,  and  estimated  the  shear  strength  at 

the  plastic  limit  to  be  100  times  that  at  the  liquid  limit,  i.e., 

2 

“  170  kN/m  at  the  plastic  limit.  The  idealized  relationship 

often  used  in  critical  state  models  (Atkinson  and  Bransby,  1978, 

Wroth  and  Wood,  1978,  and  Wroth,  1979)  is  shown  in  Fig.  2.10(a). 

Using  the  critical  state  model  and  the  above  LI  versus  c  relation- 

u 

ship,  Wroth  (1979)  developed  an  idealized  relationship  between 

vertical  effective  stress  and  liquidity  index  which  is  shown  in 

Fig.  2.10(b).  Based  on  the  experimental  evidence  Wroth  estimated 

that  o  =  6.3  kPa  at  the  liquid  limit  (or  LI  =  1)  and  o  =  630  kPa 
v  v 

at  the  plastic  limit  (LI  =  0). 

2.3.3  Normalized  Virgin  Consolidation  Lines 


Normalizing  Eq.  2.7  by  e 


L 

-C 


C  1 

— —  log  0  +  — 

P  °  v  p 

L  L 


LL 


100  s 


and  Gg  is  the  specific  gravity  of  the  soil  grains. 


(2.11) 


where  ,  the  void  ratio  at  liquid  limit,  is  obtained  from: 


(2.12) 


With  the  diffuse  double  layer  theory,  Nagaraj  and 

S r i ni va samu r t hy  (1983)  analytically  developed  and  supported  with 

experimental  data  that  the  half-space  distance,  d,  is  linearly 

* 

correlated  with  e/e^,  and  linearly  correlated  with  logOv.  This 


aspect  was  examined  further  using  the  e-logOv  curves  of  six  natural 


soils,  the  liquid  limit  water  contents  of  which  varied  over  a  four 
fold  range.  All  six  curves  fell  into  a  narrow  range,  and,  with  a 
correlation  coefficient  of  0.953,  the  best  fit  straight  line  was 
functionally  expressed  as: 

e/e  -  1.099  -  0.2237  log  o  (2.13) 

L»  1  U  V 

'  2 

where  is  in  kN/m  .  With  additional  data  for  eleven  soils, 

Nagaraj  and  Srinivasamurthy  (1986)  further  refined  the  above  equa¬ 
tion  (correlation  coefficient  of  0.962)  and  expressed  it  as: 

e/e  =  1.122  -  0.2343  log.,  o  (2.14) 

L  1  U  V 

This  was  further  confirmed  in  this  study,  using  the  critical 
state  model  and  the  data  obtained  from  Skempton  (1944)  for  16  clays 
of  different  geological  origin  (Table  2.5).  When  the  wide  scat¬ 
tered  consolidation  lines  were  normalized  by  their  corresponding 
e, '6,  they  fell  into  a  narrow  range  as  shown  in  Fig.  2.11.  The 
equation  of  the  best  fit  straight  line  was  found  as: 

e/e,  *=  1.084  -  0.  2  1  54  log,n  o  (2.15) 

L  1 0  v 

'  2 
where  o  i6  in  kN/m  . 

v 

In  addition,  conventional  consolidation  tests  were  performed 
on  three  different  artificially  sedimented  soils,  kaolinite,  "K50" 
containing  50%  kaolinite  and  50%  silt,  and  grundite,  an  illitic 
clay.  The  soil  slurries  were  made  to  a  water  content  of  approxi¬ 
mately  2  1/2  times  the  liquid  limit,  and  were  consolidated  in  a 


Table  2.5  Properties  of  Skempton's  (1944)  Clays 


No. 

Soil 

LL 

PI 

G 

s 

61 

C 

c 

eL 

eP 

1 

Boulder  Clay 

28 

14 

2.69 

0.76 

0.  12 

0.75 

0.38 

2 

Wealden  Clay 

39 

20 

2.73 

1.25 

0.24 

1 .06 

0.52 

3 

R.  Severn  Alluvium 

4  6 

2  1 

2.59 

1.22 

0.21 

1.19 

0.65 

4 

Kaolin 

50 

18 

2.64 

1.51 

0.23 

1 .32 

0.84 

5 

Oxford  Clay 

53 

26 

2.57 

1  .  56 

0.30 

1 .  36 

0.69 

6 

Belfast  Clay 

67 

37 

2.66 

1 . 64 

0.32 

1.78 

0.80 

7 

Ganges  Clay 

69 

4  1 

2.77 

2.06 

0.42 

1.91 

0.78 

8 

Gosport  Clay 

75 

46 

2.67 

2.03 

0.46 

2.00 

0.77 

9 

London  Blue  Clay 

77 

49 

2.71 

2.26 

0.49 

2.09 

0.76 

1  0 

London  Brown  Clay 

88 

56 

2.65 

2.44 

0.56 

2.27 

0.85 

1  1 

Argille  Plastique 

128 

97 

2.5  8 

3.44 

0.81 

3.30 

0.80 

1  2 

Bosporous 

36 

18 

2.71 

1 . 24 

0.25 

0.98 

0.49 

1  3 

Boston  Blue  Cl  a . 

39 

16 

2.  78 

1.22 

0.21 

1.08 

0.64 

1  4 

Vienna  Miocene 

47 

25 

2.76 

1.47 

0.30 

1  .  30 

0.61 

1  5 

Bosporous  Blue  Clay 

58 

32 

2.85 

1.61 

0.32 

1 .65 

0. 74 

i  6 

Denmark  Marine  Clay 

127 

9  1 

2.77 

4.00 

0.91 

3.52 

1 . 00 

consolidation  chamber.  Details  of  the  slurry  preparation  method, 
consolidation  chamber  etc.  are  available  in  Sivakugan  (1987)  and 
Sivakugan  et  al.  (1988b).  The  slurry  was  consolidated  to  a  verti¬ 
cal  pressure  of  10  to  20  psi  under  Kq  conditions.  At  the  end  of 
consolidation,  a  specimen  of  63.5  mm  in  diameter  and  25.5  mm  in 
height  was  trimmed  from  the  consolidated  cake  for  the  oedometer 
test.  To  ensure  full  saturation  all  oedometer  tests  were  carried 
out  under  a  back  pressure  of  30  psi.  Since  pore  pressure  measure¬ 
ments  were  not  required,  drainage  was  allowed  from  both  top  and 
bottom  of  the  specimen.  The  index  properties  and  the  parameters 
obtained  from  the  consolidation  tests  are  given  in  Table  2.6.  The 
e/e  vs.  logo  lines  resulting  from  these  experiments  fall  within 

L  V 

the  range  of  the  other  data  (Figure  2.12).  Eq.  2.14  proposed  by 
Nagaraj  and  Srini vasamur thy  (1986)  was  essentially  the  same  as  the 
average  of  the  normalized  lines,  as  shown  in  Fig.  2.12. 

2.3.4  Estimation  of  OCR 

Wroth  (1979)  used  the  critical  state  model  shown  in  Fig.  2.10 

to  predict  OCR  in  offshore  deposits,  assuming  Cr/Cc  “  0.2.  Similar 

predictions  are  possible  with  the  e/e^  ~  l°g°v  model.  From  Eq. 

2.15,  the  "unique"  (i.e.  average)  consolidation  line  can  be  plotted 

* 

in  the  e/e^  versus  logo^  plot.  This  is  the  normally  consolidated- 
virgin  compression  line  with  the  6lope  of  Cc>  By  estimating  the 
ratio  Cr/Cc  ratio,  the  6lope  of  the  swelling  line  is  known  and  it 
can  be  drawn  through  any  point  (such  as  o^  -  1000  kPa  in  Fig. 

2.13).  Then,  for  OCR  values  of  2,  5,  10,  etc.,  lines  can  be  drawn 


Table  2.6  Properties  of  Experimental  Clays 


Soil 

LL 

PL 

G 

s 

C 

c 

C 

r 

61 

eL 

Kaolinite 

63 

36 

2.65 

0.35 

0.106 

1  .93 

1 .67 

K50 

37 

23 

2.69 

0.22 

0.022 

1.17 

1 . 00 

G  rund i t  e 

51 

29 

2.75 

0.  36 

0.060 

1.76 

1 .40 
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parallel  to  the  NC  line,  through  points  with  of  500,  200,  100 
kPa,  respectively,  on  the  swelling  line.  The  rationale  for  this  is 
provided  by  the  critical  state  model,  i.e.  Kq  consolidated  and 
other  over  consolidation  lines  are  parallel  to  the  isotropic  conso¬ 
lidation  line.  Prediction  charts  are  given  in  Figs.  2.13  and  2.14, 
where  the  OCR  contours,  which  are  straight  lines  parallel  to  the 
virgin  consolidation  line,  are  shown  for  C^/C^  of  0.1  and  0.2, 
respectively. 

When  taking  samples  from  saturated  clay  deposits,  disturbance 
should  not  cause  a  significant  change  in  water  content.  Therefore, 
a  good  prediction  of  in  situ  void  ratio  at  any  depth  is  possible. 
eL  can  be  obtained  from  liquid  limit  test  on  the  remolded  soil,  and 

can  be  estimated  if  the  soil  profile  and  the  densities  of  each 

* 

6trata  are  known.  Thus,  knowing  e/eT  and  o  at  a  given  depth,  the 
corresponding  point  can  be  located  in  Figures  2.13  or  2.14  to  pro¬ 
vide  an  estimate  of  the  OCR.  The  estimate  is  not  sensitive  to  the 
assumed  value  of  the  ratio  C^/C^,  at  least  in  the  range  C.O  -  0.2. 

2.3.5  Summa  r y 

This  review  and  evaluation  of  test  data  confirmed  that,  when 
normalized  with  respect  to  the  void  ratio  at  liquid  limit,  consoli¬ 
dation  lines  fall  in  a  narrow  band.  Based  on  this,  a  simple  pro¬ 
cedure  was  described  to  predict  OCR  of  saturated  clays  using  an 
average  normalized  virgin  consolidation  line.  Since  the  OCR  esti¬ 
mate  is  not  sensitive  to  the  Cr/Cc  ratio,  a  good  OCR  prediction  is 

possible  from  the  knowledge  of  e/e^  and  o^. 
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CHAPTER  3 

MODELLING  ANISOTROPY  OF  CLAYS  -  A  NOVEL  APPROACH 

j.i  lncroduction 

Experimental  studies  conducted  in  the  past  have  demonstrated 
the  significance  of  anisotropy  on  shear  behavior  of  natural  clays 
(e.g.  Duncan  and  Seed,  1966,  Lewin,  1973,  Mitchell,  1972,  Saada, 
1976,  Saada  and  Bianchini,  1975,  Tavenas  and  Leroueil,  1977,  Ting, 
1968).  Levin's  stress  probe  experiments  (1973)  showed  differences 
between  lateral  and  axial  strains  during  isotropic  consolidation  of 
an  initially  anisotropic  specimen  of  clay.  Anisotropy  can  also  be 
induced;  Ting  (1968)  showed  that  an  anisotropically  consolidated 
sample  reaches  an  asymptotic  isotropic  state  at  an  isotropic  stress 
about  3  times  higher  than  its  initial  vertical  stress.  Experimen¬ 
tal  data  (e.g.,  Tavenas  and  Leroueil,  1977,  Graham  et  al.,  1983) 
show  that  the  yield  loci  of  natural  deposits  of  clays  are  centered 

around  the  K  line  in  triaxial  stress  space, 
o 

These  fundamental  aspects  of  soil  anisotropy  must  be  taken 
into  account  by  geotechnical  researchers  and  engineers  if  major 
improvements  in  the  interpretation  and  use  of  in  situ  measurements 
are  to  be  made.  A  simple  model  that  captures  the  basic  facets  of 
anisotropy  and  offers  closed  form  solutions  for  strength  and  other 
parameters  at  failure  along  various  stress  paths  would  be  a  valu¬ 
able  tool  to  accomplish  this.  In  addition,  the  fundamental  issues 
involved  in  in  situ  testing  can  be  best  investigated,  if  this  model 


allows  for  sopnistication  without  loss  of  the  physical  basis  of  the 
model. 

3.2  Limitations  of  the  Current  Thinking. 

Since  the  pioneering  work  of  Roscoe  and  his  co-workers  (e.g., 
Roscoe,  Schofield  and  Thurairajah  1963,  Roscoe  and  Burland  1968), 
significant  improvements  have  been  made  in  modelling  soil  behavior. 
Yet,  none  of  the  models  existing  to  date  can  oe  used  routinely  in 
the  interpretation  of  p r e s s u r erne t e r  and  especially  self  boring 
p r e s s u r erne t e r  (SBPM)  tests.  Simple  models  such  as  the  MCC  model 
(as  discussed  in  Chapte.  ?)  suffer  from  their  inability  to  capture 
the  above  features  of  anisotropy.  More  sophisticated  soil  models 
(e.g.,  Prevost,  1978,  Dafalias  and  Herrman,  1980,  and  Kavvadas , 
1982)  require  extensive  calibration  parameters  that  would  prevent 
routine  use  or  mask  the  outcome  of  the  investigations  of  in  situ 
devices.  Some  of  these  models  (e.g.,  Prevost,  1978  and  1979)  are 
limited  to  the  deviatoric  stress  space  and  thus  cannot  be  used  in 
situations  where  pore  pressures  are  involved. 

Perhaps,  more  importantly,  a  common  shortcoming  of  all  these 
models  is  that  they  do  not  provide  a  better  understanding  of  the 
physical  basis  of  the  model,  such  as  the  work  assumption  used  in 
the  isotropic  MCC  model.  With  these  shortcomings,  one  is  left  with 
the  dissatisfaction  that  the  Important  aspects  such  as  mode  of 
failure,  possibility  of  radial  cracking,  effects  of  disturbance, 
etc.,  in  SBPM  testing  can  not  be  studied  at  a  best  rational  level. 
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A  research  program  has  been  undertaken  with  the  aim  of 
addressing  the  above  shortcomings,  i.e.  of  creating  a  model  that 
captures  the  anisotropic  features  of  clays  and  i6  applicable  to 
SBPM  interpretation.  Four  stages  are  contemplated  for  this  work: 

a.  Develop  a  simple  anisotropic  theory  built  on  a  sound  physical 
basis. 

b.  Build  upon  this  model  to  a  sophisticated  level  without  loss 
of  generality. 

c.  Use  (a)  and/or  (b)  to  interpret  SBPM  data. 

d.  Use  f b )  to  study  the  fundamental  issues  of  SBPM. 

Stages  (a)  and  (c)  are  well  underway.  An  outline  of  the  essential 
features  of  the  theory  and  its  application  to  press u_emeter  prob¬ 
lems  are  reported  in  the  following;  details  are  given  in  Thevanay- 
agam  ( 1 988 a , b) . 

3.3  New  Theoretical  Developments 

The  purpose  of  this  section  is  to  provide  the  essential 
features  of  a  plasticity  model  that  captures  the  behavior  of  aniso¬ 
tropic  clays  at  the  critical  state.  The  emphasis  is  upon:  i)  The 

anisotropic  nature  of  the  yield  surface;  ii)  The  effect  of  induced 
anisotropy;  iii)  The  prediction  of  the  failure  parameters  of  clays, 
using  simple  parameters  (  4>c  and  OCR);  and  iv)  The  ability  to  inter¬ 


pret  any  measured  in  situ  strength,  and  transform  it  to  any  other 
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stress  path  applicable  to  a  given  design  problem. 

Test  results  published  in  the  literature  are  used  for  valida¬ 
tion.  These  includes  many  natural  6oils  that  have  been  tested  in 
triaxial  comp r e s s i on / e x t e ns i on  (TC/TE),  plane  strain 

comp r e s s i on /e x t e ns i on  (PSC/PSE),  model  pr e s s u re  me t e r  (PM)  and  sim¬ 
ple  shear  (SS). 

3.3.1  Background  Material:  The  Axi-symmet r i c  Model 

The  notations  are  that  of  critical  state  soil  mechanics  (Scho¬ 
field  and  Wroth  1968);  they  are  appended  to  this  chapter  (Section 
3.7).  Elastic  isotropy  is  assumed,  and  elastic  shear  strains  are 
assumed  to  exist.  The  parameters  <  and  A  refer  to  the  slope  of  the 
swelling  and  virgin  Kq  lines  in  e  vs.  Inp  plane.  All  stresses  are 
effective  stresses.  Compressive  stresses  and  strains  are  taken 
positive  and  tensor  notation  is  used.  Due  to  the  mode  of  deposi¬ 
tion,  natural  clays  are  generally  cross  anisotropic;  the  vertical 
axis  is  assumed  to  be  the  axis  of  cross  anisotropy. 

The  original  Cam  Clay  theory  (Roscoe,  Schofield  and  Thurai ra¬ 
jah,  1963)  began  with  the  work  assumption: 

p  ep  +  q  ep  «  M  p  ep  (1) 

v  q  q 

Considering  the  continuity  requirements  of  work  relation,  Eq.l  was 
modified  to  the  form  of  the  modified  cam  clay  model,  MCC,  (Roscoe 


and  Burland,  1968): 


P£v+q£n  “  P 

(ep)2  +  ( Me  p ) 2 

»  - 

V  Q 

v  q 

which,  assuming  the  associative  flow  rule,  resulted  in  an  isotropic 
yield  surface  of  the  form: 


P  "  PP, 


+  =  0 


Dafalias  (1987)  introduced  a  modification  to  the  work  assumption: 

1/2 

pep  +  qep  =  p  (ep)2  +  (M  ep)2  -r  2  aepep  ( < 

V  v  q 

where  a  is  a  no nd i me  ns i ona 1  parameter  with  an  absolute  value  less 
than  that  of  M.  Dafalias  (1987)  considered  the  term  a  to  account 

for  anisotropic  internal  residual  stresses  and  coupling  of  ep  and 

• 

ep.  However,  it  can  be  shown  that  irrespective  of  the  nature  of 
the  soil  (whether  it  is  isotropic  or  not),  Eq.4  is  strictly  true 
and  ha6  a  definite  meaning  if  the  following  relationship  is  used 
for  a  ( The va nay agam ,  1988a): 


2  a  «  2  r  +  (r2-M2)-S- 

cP 


The  selection  of  the  parameter  a  dictates  the  importance  of  the 
coupling  term  in  Eq.  4,  and  also  the  shape  of  the  yield  locus  as 
deduced  from  the  work  equation.  If  the  normality  condition  is 
assumed,  then  a  is  related  to  the  shape  of  the  yield  locus  ,  f,  by 
(from  Eq .  5 ) : 


'  •.  _  *„  %  *,  *,  ,  ■  «  %  *-  u.  *«  *-  %  J>  . 1 %  .V  %  „  •  "•  „%  _  w 


2  a  .  2  r  +  (r  -  M  )  (6) 

Thus,  in  general,  a  must  be  a  function  of  M,  p,  q  and  r  ,  i.e. 

n 

Q  ■  g1  ( M , p , q , r ^ ) .  The  assumptions  a  «  g  (r^)  (Dafalias,  1  98  7  )  or 
a  “  0  (Roscoe  and  Burland,  1968)  are  obviously  special  cases  of 
this  more  general  formulation  for  a. 

When  assuming  a  constant  value  of  a,  it  is  appropriate  to 
relate  it  to  the  end  of  Kq  virgin  consolidation  of  clay  from  a 
slurry.  At  that  stage ,  the  following  conditions  must  be  satisfied: 


e  -  e 


=  3/2  ;  —  ■= 


£  -  £ 


and  r 


4 

P 


3  (1  -  K  ) 
o 


( 1  +  2  K  ) 
o 


(7) 


where  r ^  refers  to  the  slope  of  the  virgin  Kq  line  in  p-q  plane. 


Using  the  assumption  of  elastic  isotropy  and  e  vs.  lnp  relation¬ 
ships  for  the  elastic  strains,  the  following  relationships  exist: 


P  =  (  1  +  e  )  p  .  *  e 


*//,rx  G  _  3  (1  -  2  V) 

q/(3G)  ’  K  2  (1  +  v) 


(8) 


where  K,  G  and  v  are  the  elastic  bulk  modulus,  shear  modulus  and 
Poisson'6  ratio  of  the  clay,  respectively.  Using  Eqs.  5  and  7-8,  a 


at  the  end  of  Kq  virgin  consolidation,  a^,  is  given  by 


2  2 

3  a,  r  +n  -  M 
1  'o  'o  c 


3  a 


;  wi t  h  a 


1  -  k/A 
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1  - 


K  T)o  (  l+v) 


3  A  (1  -  2v> 
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If  the  elastic  shear  strain  is  neglected  in  Eq .  7  ,  a  *l-tc/A,  and  a 

1  o 


reduces  to 


3(l-<e/A)r,  +q2-M2 
o  o  c 

3(1-k/A) 


If  the  clay  is  isotropically  consolidated  ( £  “ep*=0),  a  becomes 

q  q  o 

zero.  It  may  be  noted  that  no  assumption  on  the  shape  of  the  yield 
surface  or  on  the  flow  rule  is  necessary  to  derive  a  ^ . 

Assuming  a  constant  value  of  a,  i.e.  a  *=  a^  ,  the  shape  of  the 
yield  locus  function  is  deduced  from  the  work  equation  as: 


_J,  _  "2  ^  "'I  -  2“  ^ 

P i  M2+r2-2ar 


(10a) 


where  p  and  pi  correspond  to  any  two  positions  on  the  yield  sur¬ 
face,  f  *=  0  (Fig.  la).  However,  M  takes  the  value  (compression 

side)  for  points  located  above  point  A  in  Figure  la,  and  M  (exten¬ 
sion  side)  for  points  below  A.  Point  A  corresponds  to  the  condi¬ 
tion  of  zero  increment  in  plastic  shear  strain,  i.e.  9  f / 3  q  c  0.  At 
that  point,  the  value  of  r  is  equal  to  a  (from  Eq.  6  with  9  f  /  9  q  = 
0).  The  equations  for  the  state  boundary  surface  can  be  derived 
from  Eq.  10a  as: 

(T]_a)2  „  _  (M-a)2  +  -(-M  "Q  (q  -a)2  +  (M2-a2)  exp  -j— -  (10b) 

(M-a  ) 
c 


(-1-a-E)2-  -  (M-a)2  -E  . 

PePe  Pe  (M-a) 

c 


X-2< 

2  +  r  M2-a2 1  I  -E 


(q  -a)  +( M-a  )  - 

o  c  p 


(10c) 


parameters  pg  ,  Nq,  and  v^  are  defined  as  shown  in  Fig.  lb  and 
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Equations  10b  and  10c  are  plotted  in  Fig6.  2  and  3  for  kaolin 

using  the  following  parameters:  N  *3.23,  A*  0.19,  tc-0.06,  M  »0.85, 

o  c 

Me“-0.85,  and  r^-0.43.  The  state  boundary  surface  predicted  from 
these  equations  are  compared  to  experimental  data  on  stress  paths 
obtained  by  Atkinson  et  al.  (1987).  This  comparison  shows  that  the 
model  with  a  constant  a  captures  the  main  features  of  the  state 
boundary  surface  with  some  deviations.  The  first  deviation  occurs 
in  the  initial  part  of  the  triaxial  extension  (TE)  stress  path,  as 
this  part  of  the  TE  stress  path  lies  inside  the  state  boundary  sur¬ 
face;  however,  this  is  not  completely  unexpected,  at  least  in  the 
initial  part  of  the  stress  path,  as  the  initial  part  of  the  TE 
stress  path  remains  within  the  yield  locus  and  thus  must  be  elastic 
( The  vanay  again  and  Prapaharan,  1988).  Therefore,  this  initial  part 
of  the  TE  stress  path  is  expected  to  lie  inside  the  state  boundary 
surface.  The  stress  path  in  triaxial  compression  (TC)  closely  fol¬ 
lows  the  state  boundary  surface  up  to  the  peak  deviatoric  stress 
and  then  the  most  significant  deviation  occurs  beyond  this  point 
(clearly  illustrated  in  Figs.  2  and  3).  These  deviations  at  large 
strains  are  attributed  to  the  effects  of  induced  anisotropy. 

The  comparison  between  predictions  and  experimental  data  in 
Figs.  2  and  3  demonstrates  that  the  proposed  model  provides  a  good 
representation  of  the  material  behavior;  however,  if  the  post  peak 
behavior  in  TC  and  large  strain  behavior  in  TE  are  of  interest  then 
the  induced  anisotropy  should  be  accounted  for  in  the  model.  In 


what  follows,  the  effects  of  induced  anisotropy  are  introduced  by 
considering  that  a  does  not  remain  constant,  but  varies  during 
shearing  proportionally  to  the  initial  value  qq.  The  other  assump¬ 
tions  of  the  model,  related  to  the  other  evolution  characteristics 
of  the  yield  surface  f,  are:  (i)  the  projection  of  the  yield  sur¬ 

face  on  the  e  vs.  lnp  plane  has  a  constant  6lope  <,  and  (ii)  the 
hardening  is  such  that  point  E  (i.e.  p  in  Fig. la)  projected  on  the 
e  vs.  lnp  plane  remains  on  the  K^NC  line,  independent  of  plastic 
loading. 

The  expressions  for  undrained  strength  ratio  in  TC  and  TE  are 
developed  first.  Assuming  that  the  value  of  a  is  given  by  a  =  Yao 
(i.e.  proportional  to  )  at  failure,  the  values  of  Y  in  TC  and  TE 
are  obtained  by  calibrating  the  theoretical  expressions  for 
strength  ratio  using  experimental  data.  This  is  followed  by  the 
development  of  a  3-D  model.  The  results  of  the  axisymmetric  and 
3-D  models  are  used  to  obtain  closed  form  solutions  for  failure 
parameters  of  clays.  Predictability  is  demonstrated  by  comparison 
with  experimental  data.  Subsequently  the  application  of  this  model 
to  the  pressureme te r  will  be  illustrated. 

3.3.2  Undrained  Shear  in  Triaxial  Compression 

Referring  to  Table  1,  the  undrained  strength  in  TC  is  given 


OC  VO  ^  'OC  (TT2K 


K 

X  OCR  A;  A  -  1  -  j  (11a) 

2  2 

M  +r-2ar 
c _ o _ _£ 

2M2  -  2a  M 
c  c 

a  is  the  initial  vertical  stress  before  shearing,  a  represents 

the  anisotropic  state  at  failure,  and  nc  and  oc  stands  for  normally 

consolidated  clay  and  overconsolidated  clay,  respectively.  The 

shear  strength  is  as  defined  in  Table  1.  If  a  is  set  to  zero, 

Eq.llb  reduces  to  the  strength  expressions  for  CK^UC  strength  and 

CIUC  (when  r  -=0,f or  K  =1)  strength  of  the  MCC  model  (Wroth  1984  , 
o  o 

also  see  Chapter  2  of  this  report).  If  desired,  simplifications  to 
Eq.ll  can  be  made  as  follows. 

The  relationship  between  the  slope  of  the  virgin  consolidation 
line  and  swelling  line  in  e  vs  Inp  plane  (denoted  by  A  and  k)  and 
those  in  e  vs  lno^  plane  (denoted  by  A^  and  ky)  are  given  by  (The- 
vanayagam,  1988a): 
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Using  typical  Kq-OCR  relationships  (e.g.  Brooker  and  Ireland,  1965) 


for  clays,  Eq.12  reduces  to: 
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(13) 


k  »  1.3 


Thus,  A  in  Eq.  11. b  can  be  written  as: 
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Furthermore,  substituting  Eq.13  in  Eq.ll,  ( — - — can  be  expressed 
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For  typical  values  of  soil  parameters,  k^/i  Z  0.1  -  0.3  and  typical 
K^-OCR  relationships  (e.g.  Brooker  and  Ireland,  1965),  is  close 
to  1.0,  reducing  Eq.  15  to: 
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(16) 


The  form  of  Eq.16  is  cimilar  to  the  empirical  relationship  proposed 
by  Ladd,  et  al.  (1977)  and  the  theoretical  relationship  for  isotro¬ 
pic  clays  (Wroth  1984). 
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3.3.3  Undrained  Shear  in  Tiiaxial  Extension 


Similarly,  the  undrained  shear  strength  in  TE  is  given  by 


(see  Table  1): 


T,  M  M2-a2  Tr  M  (M  -a) 1(1- 

/  f  s  _ee  ,  f  c  c  | 

' a  ;TE  "  M  *  2  2  o  ;TC  (M  -a)  M 

vo  c  M  -a  vo  e  e 
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1  .  3< 


where  (t  /o  )  is  given  by  Eqs.ll  or  16;  in  this  case,  a 

L  VO  1  L 


represents  the  anisotropic  state  at  the  critical  state  after  shear¬ 
ing  in  TE.  Note  that  the  value  of  Mg  is  negative.  Classically, 
critical  state  parameters  and  are  obtained  using  the  expres¬ 
sions  (Roscoe  and  Burland  1968  ,  Dafalias  1987  ):  Mc*=6  s  i  n  <f>  /  3-s  i  n  <p 

and  =-6  s  i  n  p/ 3  +  s  in  <f> ,  where  the  p  angle  is  the  angle  of  friction  at 
failure  obtained  from  TC  (sin  p  B ( o , -a _ ) / ( a  +o  )  at  (o  -o_)  ). 

I  J  1  J  1  3  IDS  X 

However,  this  expression  is  valid  only  for  truly  isotropic  clays. 
Existing  data  (e.g.  Saada  et  al.,  1973)  show  that  the  friction 


angle  in  TE  is  generally  greater  than  that  in  TC.  In  this  report, 

using  the  failure  criteria  developed  later  in  Eq.23b,  a  better 

relationship  was  obtained  for  p  and  p  (  demonstrated  in  Eq.23c  ), 

c  e 

and  comparison  with  experimental  data  showed  very  good  agreement. 
Using  this  information,  Eq.17  reduces  to  a  simple  expression  (abso¬ 
lute  value  of  Tj  in  TE  is  taken): 
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3.3.4  Effect  of  Induced  Anisotropy  and  Variations  in  Soil 
Parameters 

The  effects  of  variability  on  soil  parameters  v  and  (k^/A)  on 
the  expressions  for  strength  ratio  in  TC  and  TE  were  studied.  Using 
the  experimental  data  on  many  clays  the  values  of  Y  in  TC  and  TE 
were  calibrated.  The  following  three  cases  were  considered: 

a.  a  «  0  (i.e.  isotropic  MCC  model); 

b.  Initial  value  of  a,  a^,  given  by  Eq.  9b  and  a  at  critical 

state  given  by  a  =  y  a  ; 

o’ 

c.  Case  (b),  but  using  Eq .  9a  for 

Case  (a)  when  compared  to  cases  (b)  and  (c)  is  used  to  study 
whether  the  initial  anisotropy  is  an  important  factor  in  the 
undrained  strength  ratio  of  clays.  Cases  (b)  and  (c)  were  chosen  to 
study  the  influence  of  Poisson's  ratio  on  the  strength  ratio. 

Using  these  cases,  the  values  of  Y  for  TC  and  TE  were  calibrated 
from  experimental  data  on  peak  strength. 

A  reasonable  variation  of  the  ratio  ky/A  in  the  range  0.1-0. 3 

showed  that  the  above  relationships  were  fairly  insensitive  to  this 

ratio,  and  most  calculations  were  thus  made  with  k  /A  ■  0.20.  For 

v 

the  TC  strength  (Eqs.11-16),  case  (c)  was  virtually  unaffected  by 
any  variation  in  v,  and  the  results  were  very  similar  to  that  of 
case  (b).  Based  on  regression  analysis,  case  (b)  was  found  to  be 
the  best,  thus  indicating  that  Eq.9b  is  appropriate  for  a^.  The 

single  most  influential  parameter  was  .  The  comparison  of  the 


node!  prediction  with  experimental  data  at  peak  strength  in  1C 


showed  that  a  is  not  a  very  sensitive  parameter  in  TC ,  and  thus,  as 
a  first  approximation,  Y  “1.0  is  a  reasonable  assumption  in  TC. 

The  predictability  of  residual  strength  of  clays  was  not  studied  at 


this  stage;  it  is  expected  that  Y  “  1.0  will  not  be  a  satisfactory 
assumption  if  the  strength  at  very  large  strains  is  of  interest. 
Also  the  rate  effects  were  not  considered  in  the  process  of  cali¬ 


bration. 


For  TE ,  the  influence  of  induced  anisotropy  was  significant. 


A  value  of  Y=  0.6  was  found  to  be  most  appropriate  for  the  data 
evaluated,  resulting  in  Fig. 4.  However,  since  the  term  (t ^  /  o  ) 
in  E  q .  1 8  is  not  very  sensitive  to  a  ,  for  simplicity  in  exp  re  •=■=•<  ng 
the  strength  ratio  in  terms  of  that  in  TC ,  it  can  be  sufficient  to 
substitute  a  “  ydq  only  in  the  first  term  of  Eq.18.  For  highly 
o ve r c on s o 1 i da t e d  clays,  since  most  of  the  stress  path  lies  inside 
the  yield  locus,  the  induced  anisotropy  may  be  expected  to  be 
smell,  i.e.  Y  may  tend  towards  1.0.  Further  simplifications  can  be 
made  to  Eq.18  as  follows;  Using  a  typical  value  of  k^/A  =  0.2,  tak¬ 
ing  a  Taylor  series  expansion,  and  neglecting  higher  order  terms 


Eq.18  reduces  to: 
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appear  applicable  to  many  clays. 

Using  Eqs.  11  to  19  and  q^  *  Mp^.  at  failure,  the  stresses  and 
shear  induced  pore  pressure  parameters  at  failure  can  be  derived 
for  both  TC  and  TE  (Table  2).  Table  1  contains  the  expression  for 
Skempton's  A-parameter  as  well. 

The  determination  of  failure  parameters  for  modes  of  failure 
other  than  TC  and  TE  is  presented  in  the  next  sections.  Two  Impor¬ 
tant  aspects  are  discussed:  (i)  development  of  the  necessary 

theoretical  relationships  required  to  obtain  failure  parameters  in 
any  mode  of  failure;  and  (ii)  rationale  for  failure  criteria  of  an 
anisotropic  material,  and  physical  explanation  for  the  forms  chosen 
in  this  development. 

3.3.5  General  3-D  Behavior  of  Clays  at  Critical  State 

Let  the  failure  criterion  for  a  clay,  given  by  f^  =  0, 
satisfy  the  critical  state  condition,  and  be  dependent  on  the  ini¬ 
tial  state  of  the  clay.  For  any  arbitrary  point  A  on  this  failure 

criterion,  the  corresponding  state  of  the  yield  surface  is  given 

y 

by  f  .  No  point  can  lie  outside  the  failure  criterion,  and  no 

A 

intersection  of  yield  and  failure  criteria  is  possib1e  except 
tangentially.  Using  the  associative  flow  rule,  the  gradient  of  f^ 
at  A  is  perpendicular  to  the  hydrostatic  stress  axis,  as  no  plas¬ 
tic  volumetric  strain  takes  place  at  critical  state.  Since  f^  and 

f^  are  in  contact  tangentially  at  A,  it  follows  that  the  gradient 

f  y  f 

of  f  “0  is  parallel  to  the  gradient  of  f  at  A,  where  f  •*  0  is  the 


OA 
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deviatoric  component  of  the  failure  criterion,  obtained  by  pro¬ 
jecting  ff  at  a  given  initial  state  onto  the  deviatoric  plane  in 
stress  space.  Since  f^  is  dependent  on  the  initial  state  of  the 
clay,  the  function  f^  is  also  dependent  on  the  initial  state.  The 
general  form  of  f^  in  3-D  stress  space  may  be  obtained  by  the  evo¬ 
lution  of  the  intersection  of  f^  and  another,  non  deviatoric  sur¬ 
face,  fp*=0.  In  this  development  f^  is  assumed  to  be  independent  of 
the  initial  state  of  the  clay.  The  resulting  failure  surface  f^  is 
given  by : 

ff*=  ff=  0,  with,  f^=  0,  and  ff  **  0  at  failure.  (21) 

d  p  d  p 

The  strain  increment  upon  loading  at  failure  (i.e.  at  critical 
state)  is  given  by: 


ij  ®  °i j 

where  3  is  a  proportionality  constant.  In  undrained  shearing  only 
5  of  the  6  components  of  strain  are  independent.  Therefore,  Eq.  22 
yields  only  five  independent  equations.  The  functions  f ^  »  0  and 
fp  ■  0  offer  two  more  independent  equations.  The  unknowns  are  the 
failure  stresses  (6  components)  and  the  proportionality  constant  3. 

Therefore  with  the  knowledge  of  the  strain  increment  at  failure 

• 

j  ,  and  Eqs.  21  and  22,  the  failure  stresses  o  can  be  solved  for 
in  closed  form,  i.e.  the  strain  path  and  Eqs.  21  and  22  lead  to  the 


solution  for  failure  stresses. 


The  choices  made  for  and  f  and  their  physical  meanings  are 
as  follows.  The  form  of  f^  is  given  by  a  general  Mise6  criteria 
type  : 


£r  -  k2-° 


(23a) 


with  S  =  a  -  -r  o.  .  6  ,  • 

ij  ij  3  kk  ij 


where  6,  .  is  the  Kronecker  delta,  oj  .  the  stress  tensor,  S,  .  the 
ij  ij  ’  ij 

deviatoric  stress  tensor,  and  represents  the  anisotropic  state 

of  the  material  at  failure,  with  k  and  ot^  dependent  on  the  initial 
state  of  the  material  only.  If  the  clay  is  initially  cross- 
anisotropic  (e.g.  natural  1-D  consolidated  clays)  the  axes  of 
stresses  and  strains  are  chosen  to  coincide  with  the  axes  of  aniso¬ 
tropy  (Fig. 5).  If  a  different  coordinate  system  is  chosen,  the 
results  of  this  study  should  be  used  with  appropriate  direction 
cosines  (as  later  indicated  by  Eq.  25).  For  cross  anisotropic 


clays,  a33  *  a 


22 


-  aij/2,  and  a^.-O  for  i*j. 


The  form  for  f^  is  given  by  an  extended  von  Mises  type: 


f  )  ?  7  7  u  ^  7 

fp=  (Oj-c^)  +(o2-o3)  +(o3~°i)  -—  (o1+o2  +  o3)  =  0 


(23b) 


where  °j»°2  anc*  °3  are  principal  stresses  and  H  is  a  constant.  The 
failure  criteria  given  in  Eqs.  23a  and  23b  are  interpretable  within 
the  context  of  mechanics.  The  terms  can  be  considered  to 

represent  the  deviatoric  component  of  residual  stresses  in  the 


clay,  that  is  locked-in  stresses  indicative  of  the  anisotropic 
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nature  of  the  soil.  The  deviatoric  stresses  internally  available 
in  the  soil  are  given  by  (S^  -  a^).  From  an  energy  standpoint, 
Eq.  23a  can  be  interpreted  in  terms  of  the  shear  strain  energy 
internally  stored  in  the  soil,  while  Eq.  23b  shows  that  the  strain 
energy  externally  available  to  do  work  is  dependent  on  the  mean 
stress.  Taken  together,  those  equations  imply  that  the  mean  and 
deviatoric  stresses  at  failure  are  dependent  on  the  orientation  of 
the  failure  plane. 


The  general  3-D  failure  criteria,  Eqs.  23a-b,  should  satisfy 

some  specific  conditions:  (i)  q^  *  M  p^  in  TC  and  TE ;  and  (ii)  the 

results  obtained  using  the  axisymmetric  model  given  earlier.  These 

conditions  can  be  used  to  solve  for  the  model  parameters  a  ,  k, 

2  2 

and  H.  The  first  condition  implies  that  H  »  M  ,  and  it  can  be 
used  to  obtain  the  relationship  between  and  $  : 


s  i  n  4> 


sin 
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(23c) 


The  predictions  made  with  Eq.  23c  are  given  in  Fig.  6  (dashed 
line)  for  the  experimental  data  used  earlier.  This  agreement 
between  predicted  and  experimental  values  is  very  good.  For  com¬ 
parison,  the  prediction  made  using  the  work  by  Lade  and  Mussante 
(1977  and  1978),  and  Matsuoka  and  Nakai  (1982)  are  also  given  in 
Fig.  6.  For  a  given  ,  these  relationships  tend  to  underpredict 
the  value  of  <t>e.  The  results  in  Fig.  6  indicate  that,  based  upon 


the  TC  and  TE  failure  modes,  the  proposed  criteria  are  satisf ac- 


tory.  However,  alternate  forms  could  be  explored  based  upon  data 
for  other  modes  of  failure. 


Using  the  second  condition  that  the  failure  criteria  must 
satisfy,  i.e.  the  axisymmetric  results  in  TC  and  TE ,  the  two 
unknowns  and  k  in  Eq.23a  can  be  solved  for  (Thevanayagam , 

1  988a): 
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where  (t^/o  )^£  is  given  by  Eqs.ll 

3.4  Summary  of  the  New  Model. 
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As  indicated  earlier,  with  the  knowledge  of  strain  increments 

at  failure,  Eqs.  22  to  23a-b  gives  seven  independent  equations. 

6  sin  4> 

There  are  7  unknowns  (o. ,'s  and  8).  Using  OCR,  M  =  — - - - , 

ii  c  3  -  sin  $ 

c 

Eqs.  9b,  11  or  16,  24,  a  =  y  a  ,  and  o  ,  all  the  model  parameters 

o  vo 

can  be  determined.  Consequently,  using  Eq .  22  and  23a-b,  the 
failure  stresses  in  any  mode  of  failure  can  be  solved  for  in  closed 
form.  As  noted  earlier,  if  a  coordinate  system  different  from  that 
shown  in  Fig.  5  is  chosen,  appropriate  direction  cosines  must  be 
applied  to  Eqs.  22  and  23a-b.  For  example,  if  the  strain  incre¬ 
ments  at  failure  are  known  in  a  new  coordinate  system  different 


from  that  in  Fig.  5,  Eq.  22  should  be  modified  as  follows: 
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where  o  are  stresses  in  the  coordinate  system  shown  in  Fig.  5, 

1  are  the  direction  cosines,  i.e.  cosine  of  the  angle  between 

axis  i  in  the  system  shown  in  Fig.  5  and  axis  j  in  the  new  system, 

•  * 

and  e  .  .  is  the  strain  increment  at  failure  in  the  new  coordinate 
1 J 

system. 


Using  Eqs.  15,  and  21-24,  the  strength  relationships  for  some 
of  the  most  commonly  encountered  modes  of  failure  are  given  by  (see 
Table  2): 
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with  a 


Y  a 


The  factor  0.9  in  the  above  (Eq.  26b)  was  introduced  to  take  into 
account  the  non-uniformity  in  stress  distribution  in  simple  shear 
(SS)  tests  conducted  in  the  laboratory,  which  was  found  to  reduce 
the  strength  of  the  clay  by  about  10%  (Prevost  and  Hoeg  1976). 

Additional  information  on  stresses,  friction  angles,  and  shear 

induced  pore  pressure  parameters  at  failure  are  given  in  Table  2 

for  various  stress  paths.  Eqs.  21-23  can  be  used  to  generate  the 

same  information  for  any  other  mode  of  failure  that  is  of  interest. 

Invoking  Jaky's  relationship  (1948)  for  K  (K  -■  1  -  sine*)  ),  the 

oo  c 

relationships  for  friction  angle  and  undrained  shear  strength  in 
various  modes  of  failure  can  be  further  simplified.  Figs.  7  and  8 
show  these  simplified  relationships. 


From  simple  shear  cest  data  on  natural  clays,  Ladd  et  a  1 . 
(1977)  showed  that  the  strength  ratio  of  overconsolidated  clays  can 
be  related  to  that  in  NC  clays  and  OCR  by: 


Tf  Tf 

( — )  =  ( -) 

o  'oc  o  nc 

vo  v  o 


OCR 


(29) 


where  m  is  a  constant,  approximately  equal  to  0.80.  For  isotropic 
clays,  Wroth  (1984)  showed  a  similar  theoretical  expression  for 
strength  ratio  in  TC.  Using  the  theoretical  development  reported 
herein,  it  is  possible  to  show  a  general  theoretical  relationship 
for  many  modes  of  failure  for  1-D  consolidated  clays.  Recalling 
Eqs.  1  6  ,  20  ,  26-28  ,  the  following  holds  true  for  any  particular  type 

of  undrained  shear  in  TC/TE,  PSC/PSE,  PM  and  SS: 


a  y~  (S)Cl 


(DEG) 
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VO 


A 


OCR 


(30) 


Figs.  4,  9a-b  and  10  show  the  measured  and  predicted  strength 
of  several  natural  clays,  for  many  different  stress  paths,  geologic 
origins  and  OCR.  The  close  agreement  observed  in  these  figures  as 
well  as  in  Fig.  6  substantiates  the  proposed  theory  and  the  choice 
of  the  failure  criteria.  The  general  concept  of  the  approach  and 
its  advantages  are  summarized  as  follows: 


u 


i.  The  notion  that  mean  stress  at  failure  is  dependent  on  the 
mode  of  undrained  failure  is  included; 

ii.  One  axisymmetric  yield  surface  and  a  general  form  for  the 

failure  surface  ff  are  introduced,  which  allow  the  complete 
determination  of  failure  parameters  for  any  mode  of  failure.  4 

iii.  The  initial  anisotropy  is  properly  taken  into  account  in  the 

yield  surface  and  failure  criteria.  i 


kV- 


i v .  The  effect  of  induced  anisotropy  is  captured  and  incorporated 
by  a  calibration  procedure.  Given  the  initial  conditions  of 
the  soil,  the  complete  failure  behavior  at  critical  state 
(strength,  pore  pressures,  friction  angles,  etc.)  along  any 
mode  of  failure  is  given  by  the  relationships  in  Table  2  or 
Eq  s .  2  1-23. 


v.  Limited  laboratory  testing  of  soil  is  required  to  calibrate 
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vi.  Any  relationship  derived  in  this  work  can  be  related  using 

the  two  most  commonly  used  soil  parameters,  <J>c  and  OCR  (with 

a  reasonable  value  of  <  /A  such  as  0.2). 

v 


vii.  The  model  implementation  is  simple  and  could  be  used  during 
field  tests  as  experiments  are  being  performed. 


To  the  authors'  knowledge,  it  is  the  first  time  that  a  model 


which  captures  the  essential  features  of  Initial  and  induced  aniso¬ 
tropy  ,  yet  yielding  closed  form  solutions  for  all  the  failure 
parameters  of  clays,  is  presented.  Each  of  the  steps  involved  in 
the  development  is  physically  explained  and  complemented  using 
experimental  data.  The  model  strictly  adheres  to  all  the  facets  of 
anisotropy  known  to  date.  Sophistication  of  the  model  thereof  is 
possible  and  is  being  developed. 


The  proposed  theory  can  certainly  be  applied  to  a  number  of 


geotechnical  engineering  problems  (The vanayagam ,  1988a  and  1  988b), 

such  as: 


i.  stability  analysis  of  embankments,  retaining  wa 1 1 s ,  vertical 
cuts,  etc. 


ii.  interpretation  of  vane  shear  test 


iii.  study  of  SBPM,  etc., 


The  potential  that  the  theory  has  for  SBPM  interpretation  is  intro¬ 


duced  in  the  next  section 
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Table  3:  Some  Selected  Forms  of  Failure  Criteria 


■ 

Case  I 

Case  II 

Case  III 

ff 

ff+f^  =  0 

P  d 

ff  +  f*  - 
P  d 

0 

f  f  +  f  5-0 

P  d 

ff 

p 

1/2  [ ( o 1  -  o2)^  +  (o2  -  o  ^ )  2 

Max  of  ( o 

•  °  2  ) 

Max  of  S . . 

12 

+(°1  '  °3)f)  =  M2p2 

=  M  p 

=  2/3  Mp,  i=j 

‘l 

3  /  2 ( s .  .  -  a.  .  ) 
ij  iJ 

3/2<sij  - 

Qij} 

3/ 2  (  s  .  .  -  a  .  .  ) 

12  12 

( s  .  .  -  a .  .  )-k2  =  0 

12  1J 

<su  '  “u’- 

k2  =  0 

(s  .  .  -  a .  .  )-k2  =  0 

12  il 

Since  the  induced  anisotropy  factor  Y  was  calibrated  using 
data  on  peak  strength  in  TC  and  TE ,  the  model  is  not  applicable  to 
obtain  the  residual  strength  of  clays.  Also,  calibrated  values  of  Y 
=  1.0  for  TC  and  Y  =  0.6  for  TE  were  deduced  from  experimental  data 
for  many  clays.  For  individual  clays  these  values  may  be  slightly 
different.  If  data  is  available  for  a  clay  of  interest,  then  using 
the  theoretical  development  of  this  report  the  value  of  Y  may  be 
calibrated  and  used  in  Eqs.  21-24  to  obtain  failure  parameters  in 
any  mode  of  failure  that  is  of  interest  to  the  analyst  for  that 
particular  clay. 

In  this  development  two  other  forms  were  also  considered  for 
fp.  They  are  shown  in  Table  3.  These  forms  also  give  the  relation¬ 
ship  between  <p  ^  and  <f>  g  given  by  Eq.23c.  Without  further  data  on  the 
failure  stresses  of  clays,  the  relative  merits  of  these  different 
criteria  could  not  be  studied.  Eq.23b  was  selected  because  of  its 
i n t e r p r e t a bi 1 i t y  in  terms  of  shear  strain  energy  and  mathematical 
simplicity.  Based  on  existing  data,  it  is  shown  (e.g.  see  predic¬ 
tions  in  Table  5a-b)  that  f^  in  Eq.23b  predicts  failure  stresses 
reasonably  well,  and  appears  to  perform  better  than  the  other  two 
possibilities  shown  in  Table  3.  An  experimental  program  to  study 
the  relative  merits  of  these  criteria  will  certainly  be  highly 
revealing. 

3.5  Some  Novel  Applications  to  P r e s s u r erne t e r  Problem 

Applications  of  this  theory  are  focusing  on  a  6tudy  of  the 
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following  aspects  of  self  boring  pressuremeter  measurements. 

1.  State  of  stresses  at  failure  and  mode  of  failure; 

2.  Possibility  of  radial  cracking; 

3.  Methods  to  determine  and  OCR  using  SBPM  data; 

4.  Development  of  a  solution  technique  to  interpret  disturbed 
SBPM  data; 

5.  Possible  effects  of  rate  of  testing,  creep  and  relaxation 
time  on  SBPM  data,  and 

6.  Possible  effects  of  partial  drainage. 

Promising  results  were  obtained  to  date  with  regard  to  aspect 
1,2,3  and  4.  They  are  summarized  in  the  following;  Details  can  be 
found  in  Thevanayagam  (1988b). 

3.5.1  State  of  Stress 

The  stress  path  during  the  pressuremeter  test  is  not  fully 
agreed  upon  and  much  speculation  has  been  offered  (e.g.  Prevost, 
1979,  Ladd  et  al.,  1980,  Wroth,  1984,  Ladanyi,  1977,  Wood  and 
Wroth,  1977).  A  clear  understanding  of  the  mode  of  deformation 
during  testing  is  instrumental  to  the  development  of  a  superior 
lnteroretaticn  method.  Using  the  theory  reported  herein  and  a 
"bending  mechanism",  the  mode  of  failure  in  pressuremeter  testing 
of  clays  was  identified.  The  axisymmetric  nature  of  loading  in 


pressuremeter  tests  requires  that  the  soil  elements  around  the 


probe  be  within  two  extreme  modes  of  deformation,  plane  strain  and 
plane  stress.  The  state  of  stress  and  behavior  of  clays  at  failure 
for  these  two  extreme  cases  i)  plane  strain,  =  0,  and  (ii)  plane 

stress,  =  constant,  will  be  analyzed  first.  This  information 

and  available  experimental  data  (Wood  and  Wroth  1977,  Huang  1986) 
will  then  be  used  to  develop  plausible  arguments  to  show  that  the 
assumption  of  plane  strain  =  0  is  usually  satisfactory  in  the 
interpretation  of  pressuremeter  data.  Conditions  for  which  the 
plane  strain  state  may  not  prevail  are  identified. 


For  undrained  plane  strain  shearing  with  t  =0  and  e  =-e  (sub- 

Z  j 

sequently  denoted  by  the  subscript  pm),  the  methodology  of  this 
work,  Eqs. 21-23,  yields  the  following  relationships  for  stresses 
at  failure  (also  see  Table  2): 


o-ofl  i  2  M  +0.5  a  T. 

(-.r.-j)  =  (_JL)  =  _ £ _  ( - L) 

2o  pm  vo  pm  .  i  —  .  .  vo  'TC 

vo  vo  \  |  3  ( M  +a )  vo 


o  -  a 
(^-^) 


M  -1.05a  t , 

— £ -  (— I) 

o  'TC 


^°vo  pm  \ | 3  (M  +a)  °vo 


o  -0 

- ) 


M  +  1.55  a  t 

— - -  (—1) 

v  o  'TC 


^°vo  pm  \ j  3  (M  +a)  ° vo 


where  oy-on,  °T=°22>  °e=°33 


(31a) 


(31b) 


(31c) 


Old) 


The  expression  for  (o  /t  )  is  given  by: 


( 3X2+2.25) 1/2 


o 

(“) 


Tf  pm 


Ole) 


with  X  = 


2  M  +  0.5  a 
c 


\  |  3  a 


and  the  relationships  for  o  ,  oQ  and  o  are: 

r  o  v 


0 

r 

0 

vo 

p  m 

\|3X  _t  2,25  ,  x  .  0-,  ■  o_  _f 

Mc  \  +  “  °vo  TC 


0 

V 

0 

vo 

pm 

\  3X  +  2.25 


M  +  a  o 
c  vo  TC 


(  3  1  g  ) 


°e 

0 

vo 

pm 

\  3X  +  2.25 


-  X  -  0.  5 


M  +  a  o 


(  3  1  h  ) 


vo  TC 


The  stresses  at  failure  can  also  be  obtained  using  the  other 
choices  for  f^  given  earlier  in  Table  3;  they  are  shown  in  Table  4. 
Only  results  obtained  with  the  failure  criterion  in  Eqs.  23a-b  are 
reported  in  the  following  discussion. 

Eq6.31a-c  imply  that,  independent  of  OCR  and  initial  magni¬ 
tudes  of  principal  stresses  in  horizontal  and  vertical  directions, 
the  following  inequalities  hold  true: 

o  >  o  >  a,  at  failure.  (32) 

r  v  o 

and  that  the  r-6  plane  is  always  the  failure  plane.  Figs.  11a  and 


■/."/  V.-'V. 


■  .*  >'• v-  v'  v  %■  v 


Max  of 


u 

H 


O 

H 


for  NC  and  OC 


lib  show  the  variation  of  o  /o  ,  o  /o  and  or/o 

r  vo  v  vo  6  vo 


clays 

(OCR  =  10) 

as  a 

function 

of  <p 

c 

.  Fig.  11c 

gives 

the  variation 

of  oQ 

/  o  for  various 

vo 

values  of 

OCR 

and  <)>  ,  and 
c 

for  Y 

values  of 

0.6 

and  1 

. 0 .  Table 

5a-b  , 

Figs.  12a 

and 

12b  compa  re 

measured  values 

o  f 

failure  stresses 

for 

Spestone  kaolin 

(Wood  and 

Wroth  , 

1977)  to 

those  predicted  by  Eqs.  31f-h.  Evaluation  of  these  results  lead  to 
the  following  implications: 

• 

1.  As  long  as  the  plane  strain  condition,  £^  =  0,  is  maintained, 
the  stresses  o^_  and  o  ^  become  the  major  and  minor  principal 
stresses  respectively,  irrespective  of  the  initial  state  of 
stress; 

2.  There  exists  a  critical  OCR,  called  OCR  ,  at  which  (o  /o  ) 

c  v  vo 

at  failure  remains  equal  to  1.0.  For  most  of  the  clays,  this 

0  C  R  ^  lies  in  a  narrow  range  of  1.5-2. 5  (Fig.  13).  For  a  clay 

with  OCR  <  OCR  ,  o  decreases  sharply  and  becomes  the  inter- 

mediate  principal  stress  at  failure.  The  stress  o  ^  may 

increase  or  decrease  depending  on  the  characteristics  of  the 

clay  (for  very  weak  clays  o  decreases),  and  o Q  decreases 

r  o 

becoming  the  minor  principal  stress  at  failure  ( F i g s .  1  1  a - c ) . 

If  OCR  =  OCR  ,  then  o  increases  and  o  „  decreases,  while  o 
c  r  0  v 

remains  constant  and  becomes  the  intermediate  principal 

stress.  If  OCR  >  OCR  ,  o  and  o  increase  while  oQ 

decreases.  The  increase  in  depends  on  the  value  of  OCR 

and  the  clay.  However,  the  changes  in  o  and  o  are  such  that 

v  r 

o  ^  always  becomes  the  intermediate  principal  stress.  For  any 
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Note  b:  Case  I  corresponds  to  Eqs.  23a-b  in  Sec 


Table  5b:  Predicted  and  Measured  Critical  OCR 


Predicted 

Measured  j 

2.0 

<2.5  j 

Measured  value  was  intepreted  by  the  writers  using 
the  experimental  data  from  Wood  et  al.  (1977). 
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Figure  12b  Theoretical  Versus  Measured  Stress  La*-jo  at. 

Failure  Jn  Undralne-a  P. .  Strain  a  =  q 

z 

Shear i ng  of  Clays  (Measur'd  after  Wooa  ard  l 
1977) 


value  of  OCR,  the  failure  takes  place  in  the  r-6  plane.  Thi 
is  summarized  in  Fig.  lid.  Experimental  data  obtained  by 
Wood  and  Wroth  (1977)  (Figs.  12a-b)  support  these  trends. 


3.  The  mean  principal  stress,  in  general,  does  not  remain  con¬ 
stant.  For  lightly  OC  clays  the  mean  principal  stress 
decreases  whereas  it  increases  for  highly  OC  clays; 


4.  The  angle  is  always  greater  than  the  corresponding  angle 

P  ID 

$  obtained  in  TC  (Fig. 8). 
c 

5.  Fig.  11c  indicates  that,  for  $  >  37°,  o  ^  becomes  negative. 

Also,  the  chances  for  negative  value  of  o^  increases  with  y 
approaching  1.0.  Since  Y  is  expected  to  be  close  to  1.0  for 
highly  OC  clays,  the  likelihood  for  negative  o&  is  slightly 
higher  for  OC  clays  than  for  NC  clays. 


The  solution  for  failure  stresses  o  /o  ,  oj  o  and  e  /e  in 

rvobvo  zr 

undrained  plane  stress  shear  of  clays  (a ^  =  constant,  £  ^  >  0)  can 
be  obtained  by  solving  the  following  set  of  equations: 

(2o  -  o  -  a.  -  3a,  ) 2  +  (2o  -  o  -  o  +  1.5  a  )“ 

v  r  6  11  r  v  b  1  1  ' 


+  (2o  -o  -o  +  1.5  a  )2  -  6  k2  =  0 
o  r  v  11 


(°v  -  °e) 2  +  (ov  ■  °r)2  +  (or  "  °e)2  ■ 


+  0 
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2 


0 


(33b) 


20v  -  °r  -  oe  -  30n 
2°r  K5  °  1  1 


(33c) 


where  o^j  =  c>v  ,  o^0  =  0r>  0^3  =  0g,  =  £z  (vertical  strain 

•  • 

increment),  and  e^2  =  c  >  0.  These  equations  can  be  solved  numer- 

•  • 

ically.  Fig.  14  shows  the  variation  of  £  /e  as  a  function  of  t 

2  r  c 

for  various  values  of  OCR.  This  figure  shows  that,  in  general,  the 

vertical  strain  increment  is  compressive  for  lightly  0C  clays, 

whereas  it  is  tensile  for  highly  0C  clays.  For  each  value  of  <p  , 

c 

there  exists  an  OCR  at  which  the  incremental  vertical  strain  is 
zero  (i.e.  plane  strain).  This  value  of  OCR  is  the  value  of  OCR 

c 

deduced  earlier  using  the  results  for  plane  strain  shearing  of 
clays. 


Let  us  consider  the  application  of  radial  stresses  in  a 
cylindrical  cavity  wall  as  shown  in  Fig.  15.  For  the  purpose  of 
discussion,  the  element  is  chosen  as  a  reference  point  if  end 

effects  are  neglected,  or  any  other  element  below  A^  can  be 

selected  if  end  effects  are  assumed  to  exist.  In  deep  tests,  the 
soil  above  the  level  A^  can  be  considered  as  a  thick  plate  clamped 
at  some  distance  away  from  the  probe.  Using  the  above  descriptions 
of  clay  behavior  in  plane  stress  and  plane  strain  shearing,  the 
mode  of  failure  and  the  mechanism  causing  it  in  p r e s s u r e me t e r  tests 


can  be  deduced  from  the  following  reasoning. 


Assuming  that  plane  stress  conditions  prevail,  the  results  in 

Fig.  14  indicates  that  depending  on  the  value  of  OCR,  e  can  be 

z 

tensile  (highly  OC),  zero  (at  OCR^ ) ,  or  compressive  (lightly  OC). 

Considering  the  soil  layer  above  the  level  as  a  thick  plate,  a 

compressive  strain  increment,  e^,  requires  that  the  plate  with  a 

central  hole  deflect  downwards.  Since  the  flexural  rigidity  of  the 

plate  is  large,  a  small  amount  of  deflection  requires  a  substantial 

reduction  in  the  stress  distribution  o  under  the  plate,  i.e.  in 

v 

the  horizontal  plane  at  the  level  of  A^ .  This  contradicts  the  ori¬ 
ginal  assumption  that  plane  stress  in  vertical  direction  prevails. 
The  same  argument  but  using  a  hogging  type  of  movement  can  be  used 
for  the  case  of  highly  OC  clays.  It  should  be  remembered  that  this 
argument  is  valid  only  if  the  plate  is  very  thick.  However,  in  the 
case  c.f  OCR  =  OCR^  ,  this  bending  mechanism  is  absent  (e^  =  0),  and 
the  soil  elements  satisfy  both  plane  stress  and  plane  strain  condi¬ 
tions,  independently  of  the  depth  at  which  the  test  is  performed. 
Therefore,  the  soil  elements  around  the  probe  in  pressuremeter 
testing  at  large  depths  cannot  be  expected  to  experience  the  plane 

stress  mode  of  failure  except  when  OCR  =  OCR  . 

c 

Assuming  that  plane  strain  conditions  prevail,  for  OCR  =  OCR 

c 

o  remains  equal  to  o  ,  and  the  mode  of  failure  is  both  plane 
v  1  vo 

etrain  and  plane  stress  shearing,  as  already  discussed  above.  If 
OCR  >  0  C  R  ^  ,  the  theory  reported  herein  (see  Figs,  lla-c)  and  lim¬ 
ited  experimental  data  (Wood  and  Wroth,  1977)  imply  that 
increases.  To  achieve  this,  an  upward  movement  of  the  thick  plate 


of  soil  is  necessary.  Therefore,  the  soil  elements  around  the 

probe  would  apparently  deviate  from  the  plane  strain  condition. 

However,  considering  the  high  flexural  rigidity  of  the  plate,  this 

increase  in  is  possible  with  a  very  small  upward  deflection  of 

the  chick  plate  at  level  A^;  thus,  this  upward  deflection  is  in 

fact  negligible  and  the  element  at  level  can  be  considered  to  be 

in  plane  strain.  Similar  argument  (with  downward  deflection,  and 

decrease  in  o  )  is  plausible  for  a  case  of  OCR  <  OCR  .  Therefore, 
v  c 

it  can  be  concluded  that  for  pressuremeter  testing  of  clays  at 
relatively  large  depths,  the  soil  elements  closely  experience  a 
state  of  plane  strain  shearing.  The  "bending  mechanism"  of  the 
thick  plate  of  soil  above  the  pr e s s u r erne t e r  probe  is  responsible 
for  this.  However,  the  argument  will  not  hold  true  for  testing  at 
shallow  depths  except  when  OCR  =  OCR^.  With  further  development  of 
the  model  reported  herein,  a  complete  numerical  analysis  of  this 
probl  m  is  possible  and  is  proposed  for  further  study.  Neverthe¬ 
less,  based  on  the  plausible  arguments  stated  herein  the  following 
conclusions  can  be  drawn  without  loss  of  generality.  For  a  clay 
with  OCR  =  OCR^,  the  p r e s s u r e me t e r  data  can  be  interpreted  at  any 
depth  using  the  assumption  of  plane  strain  shearing.  For  deep 
tests  the  mode  of  failure  is  very  close  to  that  in  plane  strain  and 
the  theoretical  results  reported  earlier  are  applicable.  For  shal¬ 
low  depths,  the  mode  of  failure  is  far  from  that  in  plane  strain. 


3.5.2  Interpretation  of  4  c  and  OCR  from  SBPM  Data 

It  is  noted  that  Eq.31e  depends  only  on  the  value  of  4  c  >  while 
Eq.31a  depends  on  both  4^  and  OCR  of  the  clay.  Hence,  using  these 
relationships,  a  method  to  obtain  $  and  OCR  of  the  insitu  clay 
using  SBPM  data  can  be  developed.  A  step  by  step  example  of  this 
method  is  given  in  the  following: 

1.  Using  pressuremeter  expansion  and  pore  pressure  measurements, 

obtain  ( t  )  and  corresponding  ( 0  )  at  failure, 
f  pm  r  b  r  pm 

2.  Assume  a  value  for  4  ;  Estimate  M  =6sin4  /(3-sin4  )  and  a 

c  c  c  c  o 

(from  Eq . 9  b  )  . 

3.  Use  steps  1  and  2  to  check  that  Eq .  31e  is  satisfied.  Repeat 
steps  2  and  3  with  a  different  estimate  of  4c  until  conver¬ 
gence  is  reached. 

The  value  of  OCR  may  be  deduced  from  the  value  of  4  at  con- 

c 

vergence  and  Eq.  31a.  Using  data  from  model  pressuremeter  tests 
simulated  in  a  calibration  chamber  (Chameau  et  al.,  1987),  the 
parameters  4£  and  OCR  were  obtained  with  this  procedure  for  two 
clays.  The  predicted  values  are  compared  to  those  measured  from  TC 
tests  in  Table  6.  While  the  predictions  for  normally  consolidated 
clays  is  excellent,  the  OCR  is  slightly  under  predicted  for  highly 
OC  clays  (tests  CP23  and  CP16).  This  deviation  is  attributed  to 
the  presence  of  shear  stresses  on  the  pressuremeter  membrane  due  to 
the  consolidation  around  it  and  subsequent  swelling  to  prepare 


Prediction  of  Soil  Parameters  Using  Anisotropic  Theory 


P  r  ed  i 

c  t  e  d 

Mea s  u  red 

TEST  It 
SOIL 

<J> 

y  t  c 

OCR 

d> 

Tt  c 

OCR 

CP  6 
[  K  1  0  0  ] 

19.5 

1 . 02 

20.7 

1 . 0 

CP  8 
[K50] 

20. 2 

0.99 

23.5 

1 . 0 

CP  23 
[K50] 

20.8 

7.  5 

23.5 

10.0 

CP  16 
[K100] 

23.5 

7.3 

20.7 

10.0 

[villi 


highly  OC  clays  in  the  calibration  chamber  tests. 

3.5.3  Interpretation  of  Disturbed  SBPM  Data 

The  self  boring  technique  does  produce  changes  in  the  initial 
loading  condition  and  the  initial  position  of  the  cavity  wall. 

Wroth  and  Hughes  (1973)  reported  an  outward  radial  movement  of  the 
cavity  wall  corresponding  to  about  0.5%  of  radial  strain.  Denby 
(1978)  also  reported  similar  movements  of  the  cavity  wall  well 
before  the  beginning  of  the  expansion  test.  Later  studies  using  a 
larger  size  cutting  shoe  (Benoit  and  Clough,  1986)  also  pointed  out 
this  possibility.  A  recent  integrated  study  of  the  problems 
involved  in  the  pressuremeter  tests  (Prapaharan,  1987  and  Pra- 
paharan  et  al.,  1988),  showed  that  this  initial  movement  is  the 
most  influential  when  compared  to  any  other  problems  identified 
with  the  test. 

This  section  deals  with  developing  an  interpretation  technique 
that  may  be  used  to  improve  the  current  difficulties  in  interpret¬ 
ing  p r e s s u r e me t e r  data  with  initial  movement  (i.e.  initial  distur¬ 
bance).  An  outline  of  this  new  interpretation  method  is  given  and 
illustrated  using  model  pressureueter  test  data.  The  solution  pro¬ 
cedure  requires  the  expansion  curve  and  pore  pressure  measurements 
at  the  cavity  wall  (for  ove r c o n s o  1  i d a t e d  clays  OCR  is  also 
required) . 

A. 

Let  £  refer  to  the  strain  obtained  from  measurements  in  the 
r 

disturbed  test  (i.e.  £  *  measured  strain),  and  Ac  be  the  error. 

r  o 
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The  error  Ac  in  the  measured  strain  may  be  estimated  in  the  fol- 
o 


lowing  manner  [For  illustration  purposes,  the  figures  used  herein 
were  obtained  by  applying  this  technique  to  model  pr e s s u r e me t e r 
test  data  (Test  No.  CP6)  of  Huang,  1986]: 


1.  Assume  that  the  measured  strain  £  includes  some  unknown 

r 


error  Ac  : 

o 


e  =  C  +  A  e 
r  r  o 


(  3  A  ) 


2.  For  several  values  of  Ae^,  using  Eq.  3A  and  f ^  vs.  probe 


pressure  measurements,  obtain  vs  probe  pressure  data, 


This  data  can  be  fitted  using  a  spline  function  (i.e.  one 


spline  function  for  each  Ae^). 


3. 


For  each  spline  fit  in  step  2,  using  the  method  of  deriva¬ 
tives  (e.g.  Palmer,  1972,  Ladanyi,  1972,  Baugelin,  et  al., 


1972)  obtain  the  peak  strength  ,  the  corresponding  and 


strain  at  the  cavity  wall 


These  are  shown  in  Table  7  and 
Figs.  16  and  17  for  the  data  set  considered. 


• 

4. 

For 
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6. 
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Fig. 18  (for  the  data  considered),  the  correct  value  of  Ac 


Table  7: 


Resulting  OCR, 


for  Assumed  Values  of  Error  Ae 

o 


Error 

A  e 

o 

(%) 

Tf  pm 

(kPa) 

e 

r 

(%) 

e 

r 

(%) 

0 

r 

( kP  a ) 

a 

r 

*c 

(deg.  ) 

OCR 

T  f  pm 

0.25 

117.0 

0.5 

0.25 

185.0 

1.58 

62. 

0.96 

0.15 

91.0 

0.43 

0.28 

185.5 

2 . 04 

38. 

0.85 

0.05 

7  1.0 

0.36 

0.31 

191.0 

2.70 

25. 

0.90 

1 

j  0.0 

6  4.4 

0.50 

0.50 

213.3 

3.3  1 

19. 

1 . 06 

A 

o 

U, 

60.0 

0.7 

0.75 

241.9 

4.03 

15. 

1 .26 

-0.15 

53.3 

0.8 

0.95 

260  . 

4 . 88 

12. 

1.40 

■  /  %  A 


can  be  determined 


For  the  data  in  this  example  problem,  OCR  is  equal  to  1.0,  the 


error  is  found  to  be  about  Ae^  =  ♦0.02%  from  Fig.  18.  Using  this 


value  of  Ae^,  the  measured  strain  can  be  corrected,  and  true  pres- 


suremeter  strength  and  material  parameters  deduced.  For  this  exam¬ 


ple,  they  are  given  by: 


T  =  6 4 . 4 kP a  ,  — ■ 
f  pm  ’o 


0.236,  and 


19.5  deg. 


The  value  of  f  is  in  excellent  agreement  with  that  measure  in  TC , 


about  20.7  . 


3.5.4  Existence  of  Tensile  Stress  o, 


Ladanyi  (1977)  indicated  the  possibility  of  radial  cracking 


for  frozen  soils.  Wood  and  Wroth  (1977)  observed  negative  value  of 


interpreted  circumferential  stress  ofi  for  some  p r e s s u r eme t e r  tests 


in  clays.  This  was  thought  to  lead  to  radial  cracking.  If  radial 


cracking  is  present,  then  the  mode  of  failure  may  become  that  of 


unconfined  compression  (Ladanyi  1977).  This  is  a  serious  problem, 


because  all  the  current  interpretation  methods  and  the  necessary 


relationships  to  obtain  the  s t r e s s - s t r a i n  curve  of  the  soil  are 


based  on  the  assumptions  that  the  soil  around  the  cavity  remains 


intact  and  compatibility  of  strain  holds.  If  radial  cracking  is 


present,  then  those  interpretation  methods  are  no  longer  valid. 


This  subject  has  not  been  studied  analytically  in  dejtail  except 


some  model  p r e s s u r e me t e r  tests  (Huang,  1986  and  Chameau  et  al., 
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1987)  and  simulated  plane  strain  e  =  0  tests  on  clays  in  true 
triaxial  tests  (Wood  and  Wroth,  1977).  None  of  the  simulated  tests 
or  model  p re s s u r eme t e r  tests  had  any  evidence  of  negative  circum¬ 
ferential  stress.  The  only  theoretical  study  on  mode  of  deformation 
in  SBPM  was  done  assuming  soil  as  an  isotropic-elastic  material 
that  is  governed  by  a  Mohr-Coulomb  criterion  with  o  remaining 
constant  (Wood  and  Wroth,  1977)  (Fig.  19).  Clearly  this  is  not  a 
valid  assumption  for  natural  clays. 

The  theoretical  results  reported  herein,  Fig.  11c,  indicates 
that  radial  cracking  is  possible  only  when  $  is  greater  than  about 
37°,  with  slightly  higher  susceptibility  for  highly  overconsoli¬ 
dated  clays.  Therefore,  it  appears  that  the  likelihood  of  tensile 
stresses  during  the  expansion  of  the  probe  is  very  small,  possibly 
with  the  exception  of  frozen  soils,  however,  in  this  case,  the 
theory  may  not  be  applicable.  This  suggests  that  tensile  stresses 
that  have  been  reported  based  on  current  interpretation  techniques 
may  not  be  real,  but  only  a  reflection  of  the  interpretation 
method.  A  plausible  reason  for  these  observations  is  proposed 
next. 


It  can  be  noted  that,  in  general,  whenever  a  case  of  tensile 
a  0  is  reported,  it  is  always  combined  with  an  'interpreted' 
strength  much  higher  than  expected  for  that  particular  deposit.  It 
is  likely  that,  in  many  of  these  reported  cases,  the  interpreted 
strength  is  high  because  the  measured  expansion  curve  contains  an 
error  in  the  measurements  of  radial  strain  due  to  initial  movement 


of  the  probe.  This  can  lead  to  tensile  stresses  being  computed 
erroneously  as  shown  in  the  following. 


Let  the  true  material  s t r e s s - s t r a i n  of  the  material  be  given 


by  : 


°r-°6  = 


(35a) 


Then,  the  "true"  a,  is: 

6 


a  =  o  -q ( e  ) 
o  r  r 


(35b) 


Let  the  interpreted  s t r e s s - s t r a i n  of  the  material  be  denoted  by: 


(0,-Og)1  =  q1 (Cr ) 


(36  ) 


where  is  the  measured  strain,  and  the  superscript  "i"  refers  to 

"interpreted"  parameters.  The  interpreted  circumferential  stress 
i  s  : 


0  6  "  °rV(£r) 


(  37  ) 


If  the  interpreted  relation  q  is  greater  than  the  actual  relation 
q,  at  least  for  some  range  of  strain,  then: 


,  i 


°6  <  % 


(  38 


Thus,  if  q*  is  interpreted  high,  then  it  is  possible  that  o ^ 

V 

will  reach  negative  values  while  oQ  will  in  fact  be  positive.  This 
condition  can  happen  only  when  the  interpreted  strength  is  higher 
than  expected  and,  thus,  any  observation  related  to  negative 
stresses  are  most  probably  due  to  erroneous  interpretation.  How- 


ever,  since  negative  is  a  theoretical  possibility  (though  small) 

for  soils  with  high  $  ,  model  p r e s s u r e me t e r  tests  on  such  clays  are 
needed  for  clarification  of  this  issue. 

3.6  Conclusion 

A  general  3-D  anisotropic  model  has  been  developed  to  study 
the  behavior  of  clays  at  the  critical  state.  The  model  is  simple, 
and  holds  the  promise  of  being  very  successful,  as  shown  by  com¬ 
parisons  with  available  data.  With  the  developed  relationships  one 
can  predict  the  strength  and  other  failure  parameters  in  any  par¬ 
ticular  mode  of  failure.  In  addition,  given  the  initial  condi¬ 
tions,  one  can  transform  the  data  to  predict  the  behavior  in  a  dif¬ 
ferent  mode  of  deformation.  Using  the  theoretical  development  of 
this  report  and  the  data  from  SBPM,  the  strength,  friction  angle, 
and  other  failure  parameters,  along  many  different  modes  of  failure 
other  than  that  of  SBPM  can  be  predicted.  Also,  the  in  situ  OCR 
can  be  predicted  using  SBPM  data  and  this  theory.  A  novel  approach 
has  been  presented  to  solve  the  problems  caused  by  initial  distur¬ 
bance  in  SBPM  testing.  Possibility  of  radial  cracking  is  shown  to 
be  a  rare  case.  The  theory  holds  the  promise  of  being  developed 
further.  It  offers  a  valuable  model  to  study  the  fundamental 
aspects  of  many  insitu  testing  devices  in  a  new  perspective. 

3.7  Notations 

The  following  symbols  are  used  in  this  paper: 


rvnr:nr /  w*j  rv  W'i  w\[^  *  "X n i.-1  tv; 


CIUC 


CK  UC 
o 


K  NCL 
o 


■  Skempton's  A-parameter  at  failure; 

*=  isotropically  consolidated  undrained  shear 
test  in  TC; 

=  1-D  consolidated  consolidated  undrained  shear 
test  in  TC; 

=  critical  state  line; 

=  void  ratio; 

=  f(p,q,r  )  -  yield  surface; 

=  shear  modulus; 

=  bulk  modulus; 

=  coefficient  of  earth  pressure  at  rest; 

=  1-D  consolidated  normal  compression  line; 

=  specific  volume  intercept  of  the  normal  compression 
line  -  v  vs.  lnp  at  p=l  kPa; 

=  q/p,  stress  ratio  at  critical  state; 

=  p r e s s u r e me t e r  test; 

=  plane  strain  compression; 

=  plane  strain  extension; 

=  plane  stress  shearing; 

=  mean  effective  stress; 

=  exp  (N-v/A),  equivalent  pressure; 

=  deviatoric  effective  stress; 

=  memory  variable,  hardening  parameter; 

=  simple  shear; 

=  triaxial  compression; 

=  triaxial  extension; 

=  1+e,  specific  volume; 

=  (v  +  A  lnp),  equivalent  specific  volume; 

=  anisotropic  parameter; 

=  anisotropic  parameter  at  virgin  consolid;  d  state; 
=  anisotropic  tensor; 

=  strain  tensor; 

=  volumetric  strain; 

2 

*  y  (£  -  deviatoric  strain; 

*=  slope  of  the  swelling  line 
in  e-lnp  plane ; 

*  slope  of  the  swelling  line 

in  e-1 n ( o  )  plane; 
v 

■  slope  of  the  normal  compression  and  critical  state 
lines  in  e-lnp  or  e-lno  plane; 


V 


n 


o 


ij 


“  Poisson's  ratio; 

■  q/p,  stress  ratio; 

*  (q/p)  ,  stress  ratio  during  virgin  consolidati 

=  effective  angle  of  shearing  resistance; 

3  f  /  3  p 

=  rmt 

=  stress  tensor; 


Subscripts: 

c  =  compression; 

e  =  extension; 

nc  =  normally  consolidated  clay; 

oc  =  over  consolidated  clay; 

vo  =  initial  state  in  vertical  direction; 

z  =  vertical  direction 


Superscripts  : 


e 

P 


elastic  part  of 
plastic  part  of 
increment  of  . . . 
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NOTAT  I  or; 

S"IL 

R  E  F  E  R  E  N  C  !1 

1 

Florida  clay  -  1 

o 

Florida  clay  -  2 

;  t 

Florida  clay  -  3 

, 

,  Uvdrite  -  10-2 

i  ' 

liydrite  -  10-3 

Saada,  c-t  a  1  .  (1975) 

1  6 

Oirundite  -  2 

:  : 

At  ch a  f  a  1  a  y  clay  -  1 

>•, 

A  t  c  li  a  f  a  lay  cl  a  v  -  2 

;  C, 

A  t  r  li .)  f  a  }  a  y  c  1  a  v  -  3 

]  0 

A  l  c  1 1  a  f  a  1  a  v  c  1  a  v  -  4 

]  i 

Spestonr  kaolin 

Parry  and  N  a  d  a  r  a  s  a  li  (1974) 

■  1  2 

Boston  Blue  clay 

Ladd  (1965) 
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CHAPTER  4 


CONCLUSIONS 

This  chapter  summarizes  the  main  conclusions  drawn  from  this 
research  program. 

1.  Although  CK^UC  test  and  CIUC  test  follow  different  stress 

paths,  it  has  been  shown  that  a  good  estimate  of  the  CK  UC  strength 

o 

can  be  obtained  from  a  single  CIUC  test  on  the  same  soil.  It  is  a 

useful  relation  since  the  conventional  CIUC  tests  are  much  easier 

to  run  than  CK  UC  tests, 
o 

2.  Expressions  were  derived  for  the  normalized  shear  strength 
t^/o vo  using  the  Cam  clay  model,  the  modified  Cam  clay  model,  and 
the  extended  Cam  clay  and  modified  Cam  clay  models  with  spacing 
ratio.  The  normalized  shear  strength  is  a  function  of  the  friction 
angle  and  the  consolidation  characteristics  of  the  soil.  Both 
extended  models  with  spacing  ratio  predict  values  that  compare  well 
with  experimental  data.  Unlike  the  excellent  agreement  with  exper¬ 
imental  results  obtained  for  T  / o  ,  predictions  made  for 

f  vo 

Skempton's  A  parameter  at  failure  A^  were  marginal  for  all  models. 

Ac  is  very  sensitive  to  A,  which  is  a  function  of  C  and  C  .  This 
i  c  r 

severely  influences  A  and  thus  affects  the  predictions  of  A^.. 

3.  A  simple  procedure  is  described  to  predict  OCR  of 
saturated  clays  using  normalized  virgin  consolidation  lines.  Since 
the  OCR  estimate  is  not  sensitive  to  the  Cr/Cc  ratio,  a  good  OCR 
prediction  is  possible  from  the  knowledge  of  e/e  and  o  . 


4.  Emphasis  of  this  work  was  placed  on  anisotropy  in  clays, 
leading  to  the  development  of  a  general  3-D  anisotropic  model  to 
study  the  behavior  of  clays  at  the  critical  state.  To  the  authors' 
knowledge,  it  is  the  first  model  that  captures  the  essential 
features  of  both  initial  and  induced  anisotropy  yet  yielding  closed 
form  solutions  for  all  the  failure  parameters  of  clays.  The  main 
characteristics  of  the  model  are  as  follows: 

i.  The  notion  that  mean  stress  at  failure  is  dependent  on  mode 
of  undrained  failure  is  included; 

ii.  One  axisymmetric  yield  surface  and  a  general  form  for  the 

failure  surface  f^  are  introduced,  which  allow  the  complete 
determination  of  failure  parameters  for  any  mode  of  failure. 

iii.  The  initial  anisotropy  is  properly  taken  into  account  in  the 
yield  surface  and  failure  criteria. 

iv.  The  effect  of  induced  anisotropy  is  captured  and  incorporated 
by  a  calibration  procedure.  Given  the  initial  conditions  of 
the  soil,  the  complete  failure  behavior  at  critical  state 
(strength,  pore  pressures,  friction  angles,  etc.)  along  any 
mode  of  failure  can  be  determined. 

v.  Limited  laboratory  testing  is  required  to  calibrate  the 


model. 
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5.  The  proposed  theory  is  applicable  to  a  number  of  geotechn¬ 
ical  problems;  in  this  report,  it  was  used  to  investigate  several 
important  issues  related  to  the  determination  of  in  situ  soil  pro¬ 
perties  using  the  self-boring  p r e s s u r e me t e r  (SBPM): 

(i)  The  states  of  stresses  at  failure  in  plane  strain  and 

plane  stress  were  studied  with  the  anisotropic  model.  Using  these 
results  it  was  shown  that  the  assumption  of  the  plane  strain  mode 
of  failure  is  satisfactory  for  relatively  deep  tests.  and 

become  the  major  and  minor  principal  stresses  respectively, 
irrespective  of  the  initial  state  of  stresses.  It  was  also  shown 
that  the  vertical  stress  and  mean  stress  do  not  remain  constant  in 
general,  in  contrast  to  the  predictions  made  by  isotropic-elastic 
analysis.  Currently  available  experimental  data  support  these 
findings.  All  indications  are  that  these  findings  settle  the  ques¬ 
tions  regarding  the  mode  of  deformation  and  the  role  of  vertical 
stress  during  the  expansion  of  the  cylindrical  cavity. 

(ii)  The  possibility  of  radial  cracking  (negative  circum¬ 
ferential  stress)  was  analyzed  in  detail  and  shown  to  be  unlikely; 
According  to  this  theory,  radial  cracking  is  possible  only  for 
large  friction  angles,  4>c  >  37°. 

(iii)  A  method  to  obtain  $  and  OCR  using  the  SBPM  (ideal) 
data  was  outlined,  and  illustrated  through  a  step  by  step  solution. 
Using  the  theory  developed  in  this  work,  the  SBPM  failure  parame¬ 
ters  can  be  used  to  determine  the  same  parameters  for  many  dif¬ 
ferent  modes  of  failure  other  than  that  of  the  SBPM. 
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(iv)  A  procedure  has  been  proposed  to  evaluate  the  error  in 
strain  induced  by  initial  movement  of  the  pressuremeter  cavity.  It 
is  based  upon  the  determination  of  OCR  using  the  anisotropic  model 
and  its  independent  estimation  from  other  means.  After  evaluation 
cf  the  errji  in  strain,  actual  strength  parameters  can  be  calcu¬ 
lated. 
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